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I. INTRODUCTION: "ASYMPTOTIC FREEDOM" IN A CRYOSTAT. 



19761 ) in the late 



The term "heavy fermion " was coined by Steghch, Aarts et al (jSteghch et al\ , 
seventies to describe the electronic excitations in a new class of inter-metallic compound with an 
electronic density of states as much as 1000 times larger than cop per. Since the origin al discovery of 



heavy fermion behavior in CeAl^ by Andres, Graebner and Ott ([Andres et al 



19751 ) , a diversity of 



heavy fermion compounds, including superconductors, antiferromagnets and insulators have been 
discovered. In the last ten years, these materials have become the focus of intense interest with the 
discovery that inter-metallic antiferromagnets can be tuned through a qua ntum phase transition 



into a heavy fermion state by pressure, inagnet ic fields or chemical doping (|von Lohnevsen 



von Lohnevsen et al. 



1994 



Mathur et al 



19961 : 



19981 ). The "quantum critical point" that separates the 



heavy electron ground state from the antiferromagnet represents a kind of singularity in the material 
phase diagram that profoundly modifies the metallic properties, giving them a a pre-disposition 
towards superconductivity and other novel states of matter. 

One of the goals of modern condensed matter research is to couple magnetic and electronic 
properties to develop new classes of material behavior, such as high temperature superconductiv- 
ity or colossal magneto-resistance materials, spintronics, and the newly discovered multi-ferroic 
materials. Heavy electron materials lie at the very brink of magnetic instability, in a regime where 
quantum fluctuations of the magnetic and electronic degrees are strongly coupled As such, they 
are an important test-bed for the development of our understanding about the interaction between 
magnetic and electronic quantum fluctuations. 

Heavy fermion materials contain rare earth or actinide ions forming a matrix of localized mag- 
netic moments. The active physics of these materials results from the immersion of these magnetic 
moments in a quantum sea of mobile conduction electrons. In most rare earth metals and insula- 
tors, local moments tend to order antiferromagnetically, but in heavy electron metals, the quantum 
mechanical jiggling of the local moments induced by delocalized electrons is fierce enough to melt 
the magnetic order. 

The mechanism by which this takes place involves a remarkable piece of quantum physics called 
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(Jones, 


2007; 


Kondo. 


1962, 


1964) 



which a free magnetic ion, with a Curie magnetic susceptibiUty at high temperatures, becomes 
screened by the spins of the conduction sea, to ultimately form a spinless scattering center at low 
temperatures and low magnetic fields. (Fig. [T]a.). In the Kondo effect this screening process is 
continuous, and takes place once the magnetic field, or the temperature drops below a characteristic 
energy scale called the Kondo temperature Tk- Such "quenched" magnetic moments act as strong 
elastic scattering potentials for electrons, which gives rise an increase in resistivity produced by 
isolated magnetic ions. When the same process takes place inside a heavy electron material, it 
leads to a spin quenching at every site in the lattice, but now, the strong scattering at each site 
develops coherence, leading to a sudden drop in the resistivity at low temperatures. (Fig[T](b)). 



Heavy electron materials involve the de nse lattic e ana 



often called "Kondo lattice" compounds (jPoniach 



og of the single ion Kondo effect and are 



19771 ). In the lattice, the Kondo effect may 



be alternatively visualized as the dissolution of localized, and neutral magnetic f spins into the 
quantum conduction sea, where they become mobile excitations. Once mobile, these free spins 
acquire charge and form electrons with a radically enhanced effective mass (Fig. [2]). The net effect 
of this process, is an increase in the volume of the electronic Fermi surface, accompanied by a 
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FIG. 1 (a) In the Kondo effect, local moments are free at high temperatures and high fields, but become 
"screened" at temperatures and magnetic fields that are small compared with the "Kondo temperature" Tk 
forming resonant scattering centers for the electron fluid. The magnetic susceptibility x changes from a Curie 
law X ~ y at high temperature, but saturates at a constant paramagnetic value x ^ ^ low temperatures 
and fields. (b)The resistivity drops dramatically at low temperatures in heavy fermion materials, indicating 
the development of phase coher ence between the scattering ofi^ the lattice of screened magnetic ions. (After 
([Smith and Riseboroughlll985l )) 
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profound transformation in the electronic masses and interactions. 
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FIG. 2 (a) Single impurity Kondo effect builds a single fermionic level into the conduction sea, which gives 
rise to a resonance in the conduction electron density of states (b) Lattice Kondo effect builds a fermionic 
resonance into the conduction sea in each unit cell. The elastic scattering off this lattice of resonances leads 
to formation of a heavy electron band, of width T^- 



A classic example of such behavior is provided by the inter-metallic crystal CeCuQ. Superficially, 
this material is copper, alloyed with 14% Cerium. The Cerium Ce^^ ions in this material are 
Ce'^'^ ions in a 4/^ configuration with a localized magnetic moment with J = 5/2. Yet at low 
temperatures they lose their spin, behaving as if they were Ce'^'^ ions with delocalized f-electrons. 
The heavy electrons that develop in this material are a thousand times "heavier" than those 
in metallic copper, and move with a group velocity that is slower than sound. Unlike copper, 
which has Fermi temperature of order 10,000K, that of CbCuq is of order WK, and above this 
temperature, the heavy electrons disintegrate to reveal the underlying magnetic moments of the 
Cerium ions, which manifest themselves as a Curie law susceptibility x ~ f^- There are many 
hundreds of different variety of heavy electron material, many developing new and exotic phases 
at low temperatures. 

This chapter is intended as a perspective on the the current theoretical and experimental 



understanding of heavy electron materials. There are important links 



between the material in 



this chapter, and the proceeding article on the Kondo effect by Jones (jJonesl . 



ter on quantum criticality by Sachdev (jSachdev 



2003), the chap- 



20071 ) and the perspective on spin fluctuation 



theories of high temperature superconductivity by Norman (iNorman 



20071 ) . For completeness. 



I have included references to an extensive list of review articles spanning thirty years of dis- 
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covery, 



Hewson 



1991 



i nclud ing books on the Kondo effect and heavy fermions 



19931) 



Lee et al. 



general reviews on heavy fermio n physics (|Fulde et al. 



1986 



Ott 



Gruner and Zawadowski 



1987 



1974 



Stewart 



Varma 



Cox and Zawadowski 



1988 



1999 : 



Grewe and Steghchl . 



19841 ). early views of Kondo and mixed valence physics 



19761 ) ■ the sol ution of the Kondo impurity nio del by renor - 



malization group and the strong couplin , 



the Bethe Ansatz method (Andrei et al. 



199 



cond u ctivity (ICox and Mapld. 



Riseborough . 



200ol : 



1995; 



g expansion (Nozieres and B 



1983 



Tsvelik and Wiegmanl. 



Sigrist and Ueda 



Tsunetsugu et al 



tical response in heavy fermions (IDeGiorg 



behavior and quantum criticality (iColeman et al. 



2007 



Miranda and Dobrosavlievid . 



andin 



198^), hea vy fermion super- 



1980 



Wilsonl . 



19761 ). 



1991al ). Kondo insu 



19971) . X-ray spectroscopy ([Allen et al. 



2005 



Stewartl . 



ators (lAeppli and Fiskl . 



19861 ). op- 



1999 ) and t he latest r eview s on non-Fermi liquid 



2001 



2001 



Flouquet 



2005; 



Varma et al. 



von Lohnevsen et al. 



2OO2I ). There are inevitable 



apologies, for this article is highly selective and partly because of lack of lack of space does not 



cover dynamical mean fielc 



Georges et al. 



19961 : 



theory approaches to heavy ferm ion physics (jCox and Grewe 



Jarrelll . 



19951 : 



Vidhvadhiraia et al. 



19881 : 



20031 ) ■ nor the extensiv e literature on the 



order- parameter phenomenology of heavy fermion superconductors reviewed in (jSigrist and Uedal . 



1991a). 



A. Brief History 



Heavy electron materials represent a frontier in a journey of discovery in electronic and magnetic 
materials that spans more than 70 years. During this time, the concepts and understanding have 
undergone frequent and often dramatic revision. 



In the early 1930's de Haas et al. (|de Haas et al. 



19331 ) in Leiden, discovered a "resistance 



minimum" that develops in the resistivity of copper, gold, silver and many other metals at low 
temperatures (Fig. [3]). It took a further 30 years before the purity of metals and alloys improved 
to a point wher e the r esistance minimum co uld be linked to the presence of magnetic impurities 



Clogston et al. 



Clogston et al. 



1962 



Sarachik et al 



19641 ). Clogston, Mathias and collaborators at Bell Labs 



19621 ) found they could tune the conditions under which iron impurities in Niobium 
were magnetic, by alloying with Molybdenum. Beyond a certain concentration of Molybdenum, 
the ir on im purities become magnetic and a resistance minimum was observed to develop. 



In 



1961 



Anderson formulated the first microscop ic model for the formation o f magnetic moments 



in metals. Earlier work by Blandin and Friedel (jBlandin and Friedell . 



19581 ) had observed that 



localized d states form resonances in the electron sea. Anderson extended this idea and added a 
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FIG. 3 (a) Resistance minimum in MoxNbi^x after (jSarachik et al\ . \l9Qm } (b) Temperature dependence of 
excess resistivity produced by scattering off a magnetic ion, show i ng un iversal dependence on a single scale, 



the Kondo temperature. Original data from (IWhite and Geballe , 



1979 1 



new ingredient: the Coulomb interaction between the d-electrons, which he modeled by term 

Hi = Un^ni. (1) 

Anderson showed that local moments formed once the Coulomb interaction U became large. One 
of the unexpected consequences of this theory, is that local moments develop an antiferromagnetic 



coupling with t 



(Anderson 



l e spin density of the 



1961 



s urrounding e l ectron fluid, described by the i r iterac tion 



Coqblin and SchriefFei , 



19691 : 



Kondo 



1961 



1964 



Schrieffer and Wolfi 



19661 ^ 



Hi = Ja{0) ■ S 



(2) 



where S is the spin of the local moment and cj(0) is the spin density of the electron fluid. In Japan, 



Kondo (jKondd . 



1962 ) set out to examine the consequences of this result. He found that when he 
calculated the scattering rate ^ of electrons off a magnetic moment to one order higher than Born 
approximation, 



1 

— oc 
r 



Jp + 2{Jpf\n^ 



(3) 



where p is the density of state of electrons in the conduction sea and D is the width of the 
electron band. As the temperature is lowered, the logarithmic term grows, and the scattering rate 
and resistivity ultimately rises, connecting the resistance minimum with the antiferromagnetic 
interaction between spins and their surroundings. 

A deeper understandi ng of the logarithmic term in th i s scattering rate required the renor- 



malization group concept (lAnderson and Yuva] 



Noziered . 



1976 



Nozieres and Blandinl . 



1980 



196S, 



Wilson 



197C 



1971 



Fowler and Zawadowskii 



1971 



19761 ). The key idea here, is that the physics 



of a spin inside a metal depends on the energy scale at which it is probed. The "Kondo" effect 
is a manifestation of the phenomenon of "asymptotic freedom" that also governs quark physics. 
Like the quark, at high energies the local moments inside metals are asymptotically free, but at 
temperatures and energies below a characteristic scale the Kondo temperature, 



Tk ~ Z)e-i/(2Jp) 



(4) 



where p is the density of electronic states, they interact so strongly with the surrounding electrons 
that they become scr eened into a siriglet state, or "confined " at lo w energies, ultimately forming a 



Landau Fermi liquid (iNozieres 



19761 : 



Nozieres and Blandi 



3, 



1980|). 



Throughout the 1960s and 1970s, conventional wisdom had it that magnetism and supercon- 
ductivity are mutually exclusive. Tiny concentrations of magnetic produce a lethal suppression 



of superconductivity in conventional metals. Ear^ 



and superconductivity by Maple et aJ. ([Maple et al 



y wor k on the interplay of the Kondo effect 



19721 ). did suggest that the Kondo screening 



suppresses the pair breaking effects of magnetic moments, but the implication of these results 
was only slowly digested. Unfortunately, the belief in the mutual exclusion of local moments and 
superconductivity was so deepl y ingrained, t h at th e first observation of superconductivity in the 



"local moment" metal UBeis ([Bucher et al 



19751 ) was dismissed by its discoverers as an a rti- 



19751 by 



fact produced by stray filam ents of uranium. Heavy electron metals were discovered in 
Andres. Graebner. and Qtt . who observed that the inter-metallic CeAl^ forms a metal in which 
the Pauli susceptibility and linear specific heat capacity are about 1000 times larger than in con- 
ventional metals. Few believed their speculation that this might be a lattice version of the Kondo 
effect, giving rise in the lattice to a narrow band of "heavy" f-electrons. The discovery of supercq n- 



1976|) , 



ductivity in CeCu2Si2 in a similar f-electron fluid, a year later by Steglich (jSteglich et al. 
was met with widespread disbelief. All the measurements of the crystal structure of this material 
pointed to the fact that the Ce ions were in a Ce^^ or 4/^ configuration. Yet this meant one local 
moment per unit cell - which required an explanation of how these local moments do not destroy 
superconduc tivity, bu t rath er, are part of its formation. 



Doniach ( Doniach 



Kondo lattice (jKasuval . 



19771). made the visionary proposal that a heavy electron metal is a dense 



19561 ). in which every single local moment in the lattice undergoes the 



Kondo effect (Fig. [2]). In this theory, each spin is magr ietically screen ed by the conduction sea. 



19851 ). was whether there could 



One of the great concerns of the time, raised by Nozieres (jNozieresl . 
ever be sufficient conduction electrons in a dense Kondo lattice to screen each local moment. 
Theoretical work on this problem was initially stalled, for wont of any controlled way to compute 
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properties of the Kondo lattice. In the early 1980's, Anderson (jAndersonl . Il98ll ) proposed a way 
out of this log-jam. Taking a cue from the success of the 1/5" expansion in spin wave theory, and the 
expansion in statistical mechanics and particle physics, he note that the large magnetic spin 
degeneracy = 2j + 1 of f-moments could could be used to generate an expansion in the small pa- 
rameter 1/A about the limit where N o o. Anderson's idea prompted a renaissance of theoretica. 



development ( Auerbach and Levin 



198 



Ramakrishnan 



1981 



1986 



Coleman 



Read and Newns 



1983a 



1983 



1987a 



Gunnarsson and Schonhammei 



making it possible to compute the X-ray 



absorption spectra of these materials and, for the first time, examine how heavy f-ba nds form within 



the Kondo lattice. By the mid eigh ties, the first de Haas van Alphen experiments (iReinders et al. 



1986; 



Taillefer and Lonzarich 



19881 ) had detected cyclotron orbits of heavy electrons in CeCuQ and 



UPt^. With these developments, the heavy fermion concept was cemented. 



On a separate experimental front, in 1983 Ott, Rudiger, Fisk and Smith (lOtt et al 



1981 



1984 ) returned to the material UBei^, and by measuring a large discontinuity in the bulk specific 
heat at the resistive superconducting transition, confirmed it as a bulk heavy electron supercon- 
ductor. This provided a vital independent confirmation of Steglich's discovery of heavy electron 
superconductivity, assuaging the old doubts and igniting a huge new interest in heavy electron 
physics. The nuniber of h eavy electron metals and superconductors grew rapidly in the mid 1980s 



(jSigrist and Uedal . 



1991bl ). It became clear from specific heat, NMR and ultrasound experiments on 



heavy fermion superconductors that the gap is anisotropic, with lines of nodes strongly suggesting 
an electronic, rather than a phonon mechanism of pairing. These discoveries prompted theorists 
to return to earlier spin fluctuation-mediated models of anisotropic pairing. In the early summer 
of 1986, three new theor etical papers were re ceived by Physical Review, the first by Beal Monod, 



Bourbonnais and Emery (jMonod et al 



19861 ) worki ng in Orsay, Fra,nce, fo llowed closely (six weeks 



later) by papers from Scalapino, Loh and Hirsch ( Scalapino et al 



19861 ) at UC Santa Barbara, 



19861 ^ at Bell Labs, New Jersey. 



California, and Varma, Schmitt-Rink and Miyake (|Mivake et al. 
These papers contrasted heavy electron superconductivity with superfluid He — 3. Whereas He — 3 
is dominated by ferromagnetic interactions, which generate triplet pairing, these works showed 
that in heavy electron systems, soft antiferromagnetic spin fluctuations resulting from the vicinity 
to an antiferromagnetic instability would drive anisotropic d-wave pairing (Fig. |4]). The almost 
coincident discovery of high temperature superconductivity the very same year, 1986, meant that 
these early works on heavy electron superconductivity were destined to exert huge influence on 
the evolution of ideas about high temperature superconductivity. Both the RVB and the spin- 
fluctuation theory of d-wave pairing in the cuprates are, in my opinion, close cousins, if not direct 
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FIG. 4 Figure from (jMonod aZj . Il986r ). one of three path-breaking papers in 1986 to hnk d-wave pairing 
to antiferromagnetism. (a) is the bare interaction, (b) and (c) and (d) the paramagnon mediated interaction 
between anti-paraUel or paraUel spins. 

descendents of these early 1986 papers on heavy electron superconductivity. 

After a brief hiatus, interest in heavy electron physics re-ignited in the mid 1990's with the 
discovery of quantum critical points in these materials. High temperature superconductivity intro- 
duced many important new ideas into our conception of electron fluids, including 

• Non Fermi liquid behavior: the emergence of metallic states that can not be described as 
fluids of renormalized quasiparticles. 

• Quantum phase transitions and the notion that zero temperature quantum critical points 
might profoundly modify finite temperature properties of metal. 

Both of these effects are seen in a wide variety of heavy electron materials, providing an vital 
alternative venue for research on these still unsolved aspects of interlinked, magnetic and electronic 
behavior. 

In 1995, Hilbert von Lohneyson and collaborators discovered that by alloying small amounts 
of gold into CeCuQ that one can tune CcCuq-xAux through an antiferro magnetic quantuni crit' 



ical point, and then rever se the process by the application of pressure ()von Lohneysen 



von Lohneysen et al. 



19961 : 



1994 ). These experiments showed that a heavy electron metal develops "non- 
Fermi liquid" properties at a quantum critical point, including a linear temperature dependence of 
the resistivity and a logari thmic dependence of the specific heat coefficient on temperature. Shortly 



thereafter, Mathur et al. (iMathur et al 



19981 ). at Cambridge showed that when pressure is used 
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to drive the antiferromagnet Celn^ through through a quantum phase transition, heavy electron 
superconductivity develops in the vicinity of the quantum phase transition. Many new examples 
of heavy electron system have co me to light in the last few years which follow the same pattern. 



In one fascinating development, (jMonthoux and Lonzarich 



19991 ) suggested that if quasi-two di- 



mensional versions of the existing materials could be developed, then the superconducting pairing 
would be less frustrated, leading to a higher transition temperature. This led experimental groups 
to explore the effect of introducing layers into the material Celn^, leading to the discovery of the 



so called 1 — 1 — 5 compounds, in which an XIno 



compound. (jPetrovic et al 



2001 



Sidorov et al. 



CeCoIn^ and most recently, PuCoGa^ (jSarrao et al. 



ayer has been introduced into the original cubic 



20021) . T wo notable members of this group are 



3). 



20021 ) . The transition temperature rose from 



0.5-ftr to 2.bK in moving from Celn^ to CeCoIn^. Most remarkably, the transition temperature 
rises to above 18K in the PuCoGa^ material. This amazing rise in Tc, and its close co nnection with 



quan tum criticality, are very active areas of research, and may hold important clues (jCurro et all 



20051 ) to the ongoing quest to discover room temperature superconductivity. 



B. Key elements of Heavy Fermion Metals 

Before examining the theory of heavy electron materials, we make a brief tour of their key 
properties. Table A. shows a selective list of heavy fermion compounds 

1. Spin entropy: a driving force for new physics 



The properties of h eavy fermion compounds derive from the par t ially filled f orb i tals o f rare 



earth or actinide ions (Fulde et al 



Stewart 



19881 : 



Grewe and Steglich 



1991 



Lee et al 



19861 : 



Ott 



19871 : 



1984 ) . The large nuclear charge in these ions causes their f orbitals to collapse inside the 



inert gas core of the ion, turning them into localized magnetic moments. 
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Table. A. Selected Heavy Fermion Compounds. 



Type 


Material 




T r B 

^Cl -^ci 


Properties 


P 


In 

mJmol~^K~' 


Ref. 


Metal 


CeCuQ 


WK 


- 


Simple HF 
Metal 




1600 


[1] 


Super- 
conductors 


CeCu2Si2 


20K 


rc=0.17K 


First HFSC 


r2 


800-1250 


[2] 


UBen 


2.5K 


rc=0.86K 


Incoherent 
metal^HFSC 


Pc ~ 

150;ui7cm 


800 


[3] 


CeCoIn^ 


38K 


rc=2.3K 


Quasi 2D HFSC 


T 


750 


[4] 


Kondo 
Insulators 


CesPtiBi-i 


~ 80K 


- 


Fully Gapped 
KI 




- 


[5] 


CeNiSn 


~ 20K 


- 


Nodal KI 


Poor 
Metal 


- 


[6] 


Quantum 
Critical 


CeCue-xAux 


To ~ WK 


=0.1 


Chemically 
tuned QCP 


T 




[7] 


YbRh2Si2 


To ~ 2AK 


Bi==0.06T 
B||=0.66T 


Field-tuned 
QCP 


T 




[8] 


SC + 
other 
Order 


UPd2Ak 


llOK 


TAi.=14K, 


AFM + HFSC 


rp2 


210 


[9] 


U Ru2Si2 


75K 


ri = 17.5K, 
T,c=1.3K 


Hidden Order & 
HFSC 


J<2 


120/65 


[10] 



Unless otherwise stated, T* denotes the temperature of the maximum in resistivity. Tc, Xc and Be 
denote critical temperature, doping and field, p denotes the temperature temperature dependence in the 



normal state. 7,1 = Cv /T is the specific heat coefficient in the normal state. [1] ( 


Onuki and Komatsubara 


1987 


Stewart et al. 


1984b). [2] 


(Geibel et al 


. 1991; 


Geibel. et al.. 


1991 


; Stedich et al. 1976). [3] ( 


Ott et al. 


1983 


1984), [4] 


Petrovic et al.. 


2001; 


Sidorov et al. 


2002). [5] (Bucher et al. 


1994; 


Hundlev et al. 


1990 


),[6] 


(Izawa et al.. 


1999; 


Pakabatake et al.. 


1992, ] 


-990). [7] (von Lohnevsen, 


1996; 


von Lohnevsen et al. 


1994 




(Custers et a 




200c 


: Gegenwart et al 


.. 2005 


Paschen et al. 2004 


Trovarelli et al. 


. 2000). [9] (Geibel et al 
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-^LfT = Spin Eaitropv (T) 



i.KI 



FIG. 5 Showing the specific heat coefficient of UBei^ after (jOtt et all Il985f) . The area under the Cv/T 
curve up to a temperature T provides a measure of the amount of unquenched spin entropy at that tem- 
perature. The condensation entropy of heavy fermion superconductors is derived from the spin-rotational 
degrees of freedom of the local moments, and the large scale of the condensation entropy indicates that spins 
partake in the formation of the order parameter. 



1991; 


Sato et al.. 


2001 


Tou et al.. 


1995 


), [10] ( 


Kim et al.. 


2003; 


Palstra et al.. 


1985) 



Moreover, the large spin-orbit coupling in f-orbitals combines the spin and angular momentum of 
the f-states into a state of definite J and it is these large quantum spin degrees of freedom that lie 
at the heart of heavy fermion physics. 

Heavy fermion materials display properties which change qualitatively, depending on the tem- 
perature, so much so, that the room temperature and low temperature behavior almost resemble 
two different materials. At room temperature, high magnetic fields and high frequencies, they 
behave as local moment systems, with a Curie law susceptibility 



X = ^ M-' = {gj^jiB?J{J + l) 



(5) 



where M is the magnetic moment of an f state with total angular momentum J and the gyromag- 
netic ratio gj. However, at temperatures beneath a characteristic scale we call T* (to distinguish 
it from the single-ion Kondo temperature TV), the localized spin degrees of freedom melt into the 
conduction sea, releasing their spins as as mobile, conducting f-electrons. 

A Curie susceptibility is the hallmark of the decoupled, rotational dynamics of the f-moments. 



13 



associated with an unquenched entropy of S = kslnN per spin, where = 2 J + 1 is the spin- 
degeneracy of an isolated magnetic moment of angular momentum J. For example, in a Cerium 
heavy electron material, the 4/^ (L = 3) configuration of the Ce^"^ ion is spin-orbit coupled into a 
state of definite J = L — S = b/2 with = 6. Inside the crystal, the full rotational symmetry of 
each magnetic f-ion is often reduced by crystal fields to a quartet (A^ = 4) or a Kramer's doublet 
N = 2. At the characteristic temperature T*, as the Kondo effect develops, the spin entropy is 
rapidly lost from the material, and large quantities of heat are lost from the material. Since the 
area under the specific heat curve determines the entropy, 

S{T) = £ ^dT', (6) 

a rapid loss of spin entropy at low temperatures forces a sudden rise in the specific heat capacity. 
Fig. [5] illustrates this phenomenon with the specific heat capacity of UBei^. Notice how the specific 
heat coefficient Cy/T rises to a value of order IJ / mol / K'^ , and starts to saturate at about IK, 
indicating the formation of a Fermi liquid with a linear specific heat coefficient. Remarkably, just 
as the linear specific heat starts to develop, UBeis becomes superconducting, as indicated by the 
large specific heat anomaly. 



2. "Local" Fermi liquids with a single scale 



The standa rd theoretical framework for describing metals is Landau Fermi liquid theory 



(jLandau . 



19571 ). according to which, the excitation spectrum of a metal can be adiabatically con- 
nected to those of a non-interacting electron fluid. Heavy Fermion metals are extreme examples of 
Landau Fermi liquids which push the idea of adiabaticity into an regime where the bare electron 
interactions, on the scale of electron volts, are hundreds, even thousands of times larger than the 
millivolt Fermi energy scale of the heavy electron quasiparticles. The Landau Fermi liquid that 



develops in these materials s 



lares much in common with the Fermi 



isolated magnetic impurity (jNozieres 



1973; 



Nozieres and Blandinl . 



iquid that develops around an 
IQSol ). once it is quenched by 



the conduction sea as part of the Kondo effect. There are three key features of this Fermi liquid: 

• Single scale T* The quasiparticle density of states p* ~ 1/T* and scattering amplitudes 
^ko-k'cr' ^ T* scale approximately with a single scale T* . 



Almost incompressible. Heavy electron fluids are "almost incompressible", in the sense 
that the charge susceptibility Xc = dNe/dfi « p* is unrenormalized and typically more 
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than an order magnitude smaller than the quasiparticle density of states p* . This is because 
the lattice of spins severely modifies the quasiparticle density of states, but leaves the charge 
density of the fluid ne(/x), and its dependence on the chemical potential /x unchanged. 



Local. Quasiparticles scatter when in the vicinity of a local moment, giving rise to a sma. 



momentum dependence to the Landau scatt ering amplitudes. (jEngelbrecht and Bedelll . 



Yamadal . 



1975 



Yoshida and Yamada 



I975I I . 



19951 : 



Landau Fermi liquid theory relates the properties of a Fermi 



liquid to the density of st ates of the 



quasiparticles and a small number of interaction parameters (jBavm and Pethickl . |1992| ) If E'ko- is 
the energy of an isolated quasiparticle, then the quasiparticle density of states p* = X^^o- "^(-^ko- — m) 
determines the linear specific heat coefficient 



7 = LimT^o ( 



(7) 



In conventional metals the linear specific heat coefficient is of order 1 — 10 mJ mol ^. In a 



system with quadratic dispersion, E]^ 
m* are directly proportional 



2m 



'- , the quasiparticle density of states and effective mass 



m 



(8) 



where kp is the Fermi momentum. In heavy fermion compounds, the scale of p* varies widely, 
and specific heat coefficients in the range 100 - 1600 mJ mol-^K-^ have been observed. From 
this simplified perspective, the quasiparticle effective masses in heavy electron materials are two 
or three orders of magnitude "heavier" than in conventional metals. 

In Landau Fermi liquid theory, a change (Jn^'o-' in the quasiparticle occupancies causes a shift 
in the quasiparticle energies given by 



<J-Ekcr = ^ /k(7,kcr'5?^k'cr' 



(9) 



k'cr' 



In a simplified model with a spherical Fermi surface, the Landau interaction parameters only 
depend on the relative angle Oi^^y between the quasiparticle momenta, and are expanded in terms 
of Legendre Polynomials as 



/k.,k.' = \ + l)PK^k,k')[i^f + era' Ft]. 



(10) 



The dimensionless "Landau parameters" F^''^ parameterize the detailed quasiparticle interactions. 
The s-wave (/ = 0) Landau parameters determine the magnetic and charge susceptibility of a 
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Landau Fermi liquid are given by (jBavm and Pethick 



1992 : 



Landau 



19571) 



Xs 
Xc 



Pb 



2 P 2„* 



eV [l-Al] 



(11) 



where, the quantities 



/,s,a 
^0 



F 



(12) 



1 + Fq'" 

are the s-wave Landau sc attering amphtudes in the charge (s) and spin (a) channels, respectively 



( Bavm and Pethickl . 



1992)- 



The assumption of local scattering and incompressibility in heavy electron fluids simplifies the 
situation, for in this case only the / = components of the interaction remain and the quasiparticle 
scattering amplitudes become 



(13) 



Moreover, in local scattering the Pauli principle dictates that quasiparticles scattering at the same 
point can only scatter when in in opposite spin states, so that 



A 



(0) 

TT 



0. 



(14) 



and hence A^ = —A'^. The additional the assumption of incompressibility forces Xc/(e^P*) << 1, 
so that now Aq = —Aq ~ 1 and all that remains is a single parameter p* . 



T his line of reasoning, first dev eloped for the single impurity Kondo model by Nozieres (jNozierea . 



1976 



Nozieres and Blandin 



Engelbrecht and Bedelll . 



1980 ) , later extended to a bulk Fermi liquid by Bedell and Engelbrecht 



19951 ). enables us to understand two important scaling trends amongst 



heavy electron systems. The first consequence, deduced from (fTTl). is t hat t he dimensionless Som- 



merfeld ratio, or "Wilson ratio" W 



^ « 2. Wilson ( Wilson 



19761 ). found that this ratio 



is almost exactly equal to two in the numerical renormalization group treatment of the impurity 
Kondo model. T he connection b e tween this ratio and the loca l Fermi liquid theory was first identi- 



fied by Nozieres, (jNoziered . 



1976 



Nozieres and Blandin. 



1980 ). In real heavy electron systems, the 



effect of spin orbit coupling slightly modifies the precise numerical form for this ratio, nevertheless, 
the observation that ~ 1 over a wide range of materials in which the density of states vary by 
more than a factor of 100, is an indication of the incompressible and local character of heavy Fermi 
liquids (Fig. [6]). 
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FIG. 6 Plot of linear s pecific heat coefficient vs Pauli susceptibility to show approximate constancy of Wilson 



ratio. (After B. Jones (|Lee et aLl . Il986f )l. 



A second consequence of locality appears in the transport properties. In a Landau Fermi liquid, 
inelastic electron-electron scattering produces a quadratic temperature dependence in the resistivity 



(15) 



In conventional metals, resistivity is dominated by electron-phonon scattering, and the coeffi- 
cient is generally too small for the electron-electron contribution to the resistivity to be observed. 
In strongly interacting metals, the A coefficient becomes l arge, and in a beautiful piece of phe- 



nomenology, Kadowaki and Woods ([Kadowaki and Woods 



19861 ) , observed that the ratio of A to 



the square of the specific heat coefficient 7^ 

A 



OiKW 



(1 X 10"'')/LiJ]cm[mol K7mJ] 



(16) 



is approximately constant, over a range of A spanning four orders of magnitude. This too, can be 
simply understood from local Fermi liquid theory, where the local scattering amplitudes give rise 
to an electron mean-free path given by 



1 



constant + 



(17) 



kpl* ' (T*)2' 

The "A" coefficient in the electron resistivity that results from the second-term satisfies A oc 



7 oc 7 . A more detailed calculation is able to account for the magnitude of the Kadowaki 
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FIG. 7 Approxima te constancy of the Kadowaki Woods ratio, for a wide range of heavy electrons, after 
(jTsujii et aLl . l2005[ ). When spin-orbit effects are taken into account, the Kadowaki Wood ratio depends on 
the effective degeneracy N = 2J + 1 of the magnetic ion, which when taken into account leads to a far more 
precise collapse of the data onto a single curve. 



Woods constant, and its weak residual dependence on the spin degeneracy = 2 J + 1 of the 
magnetic ions (see Fig. [71). 

The approximate validity of the scaling relations 

X ^ 

— PS cons, ^ cons (18) 

7 T 

for a wide range of heavy electron compounds, constitutes excellent support for the Fermi liquid 
picture of heavy electrons. 

A classic signature of heavy fermion behavior is the dramatic change in transport properties 
that accompanies the development of a coherent heavy fermion band structure (Fig. [6]). At high 
temperatures heavy fermion compounds exhibit a large saturated resistivity, induced by incoherent 
spin-flip scattering of the conduction electrons off the local f-moments. This scattering grows as the 
temperature is lowered, but at the same time, it becomes increasingly elastic at low temperatures. 
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Development of coherence in in Cei-xLaxCuQ after Onuki and Komatsubara 



( Onuki and Komatsubara . 



1987 ). 



This leads to the development of phase coherence, the f-electron spins. In the case of heavy fermion 
metals, the development of coherence is marked by a rapid reduction in the resistivity, but in a 
remarkable class of heavy fermion or "Kondo insulators" , the development of coherence leads to a 
filled band with a tiny insulating gap of order Tk- In this case coherence is marked by a sudden 
exponential rise in the resistivity and Hall constant. 

The classic example of coherence is provided by metallic CbCuq, which develops "coherence" 
and a maximum in its resistivity around T = 10 K. Coherent heavy electron propagation is readily 
destroyed by substitutional impurities. In CcCuq, Ce^^ ions can be continuously substituted with 
non magnetic La?^ ions, producing a continuous cross-over from coherent Kondo lattice to single 
impurity behavior (Fig. [8] ) . 

One of the importa nt principles of th e Landau Fermi liquid is the Fermi surface counting rule. 



or Luttinger's theorem ljLuttingei 



3). 



1960). In non interacting electron band theory, the volume of 



the Ferm i surface counts the number o f conduction electrons. For interacting systems this rule 



survives (jMartinl . 



19821 : 



Oshikawa 



2000), with the unexpected corollary that the spin states of the 



screened local moments are also included in the sum 

[ne + n 



spmsj 



(19) 



(27r)3 

Remarkably, even though f-electrons are localized as magnetic moments at high temperatures, in 
the heavy Fermi liquid, they contribute to the Fermi surface volume. 

The most most direct evidence for the large heavy f- Fermi surfaces derives from de Haas 
van Alphen and Shubnikov de Haas experiments that measure the oscillatory diamagnetism or 
and resistivity produced by coherent quasiparticle orbits (Fig. [9|). These experiments provide 
a direct measure of the heavy electron mass, the Fermi surface geometry and volume. Since the 
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FIG. 9 (a) Fermi surface of UPt^ cal c ulate d from band-theor y assuming itinerant 5/ electrons 



(Norman et al 



1988 



Qguchi and Freemanl . 



1985 



Wang et al 



that are identified by dHvA measurements, after 



Kimura et al 



1987), s howing three orbits (cr, u; and r) 



19981 ). (b) Fourier transform of dHvA 



oscillations identifying a, w and r orbits shown in (a). 



pioneering measure ments on CeCuR and U Pt^ by Reinders and Sp r ingfo r d, Taillefer anc 



in the mid-eighties (jReinders et al 



1983; 



Taillefer and Lonzarich 



extensive number of such measurements have been carried out ( Julian et al 



1998 



McCollam et al 



2005 



Onuki and Komatsubaral . 



1988 



Taillefer et al 



1992 



;onzarich 



1987), an 



Kimura et al 



19871 ). Two key features are observed: 



• A Fermi surface volume which counts the f-electrons as itinerant quasiparticles. 

• Effective masses often in excess of one hundred free electron masses. Higher mass quasi- 
particle orbits, though inferred from thermodynamics, can not be observed with current 
measurement techniques. 

• Often, but not always, the Fermi surface geometry is in accord with band-theory, despite 
the huge renormalizations of the electron mass. 

Additional confirmation of the itinerant nature of the f-quasiparticles comes from the observa- 
tion of a Drude peak in the optical conductivity. At low temperatures, in the coherent regime, 
an extremely narrow Drude peak can be observed in the optical conductivity of heavy fermion 
metals. The weight under the Drude peak is a measure of the plasma frequency: the diamagnetic 



response of the heavy fermion metal. T his is foun c 



mass enhancement of the quasiparticles (iDeGiorgi 



to be extremely small, d epressed by the large 



19991 : 



Millis et al 



1987al ) 



-aqp{uj) 



ne 
m* 



(6) 



'\u^\<Tk ^ 

Both the optical and dHvA experiments indicate that the presence of f-spins depresses both the 
spin and diamagnetic response of the electron gas down to low temperatures. 
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II. LOCAL MOMENTS AND THE KONDO LATTICE 



A. Local Moment Formation 



1. The Anderson Model 



We begin with a discussion of how magnetic moments form at high temperatures, and how they 



are screened again at low tempera tures to fo r m a 



formation is the Anderson model ( Anderson 



1961 



^ermi liquid. The basic model for local moment 



resonance 



k,a 



k,a 



(20) 



^atomic 

where Hatomic describes the atomic limit of an isolated magnetic ion and Hresonance describes the 
hybridization of the localized f-electrons in the ion with the Bloch waves of the conduction sea. For 
pedagogical reasons, our discussion will initially focus on the case where the f-state is a Kramer's 
doublet. 

There are two key elements to the Anderson model: 

• Atomic limit. The atomic physics of an isolated ion with a single / state, described by the 
model 



H, 



atomic 



EfUf + Uuf^Ufi. 



(21) 



Here Ef is the energy of the / state and U is the Coulomb energy associated with two 
electrons in the same orbital. The atomic physics contains the basic mechanism for local 
moment formation, valid for f-electrons, but also seen in a variety of other contexts, such as 
transition metal atoms and quantum dots. 

The four quantum states of the atomic model are 

1/2) E{f) = 2Ef + U 

I/O) E{f) = 

\j ' yj J , ^22) 



non-magnetic 



l/M), l/U) E{f^)=Ef. magnetic. 

In a magnetic ground-state, the cost of inducing a "valence fluctuation" by removing or 
adding an electron to the /^ state is positive, i.e. 



removm, 



g: E{f) - E{f) = -Ef>0 ^ U/2 > Ej + U/2, 
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FIG. 10 Phase diagram for Anderson Impurity Model in the Atomic Limit. 

adding: E{f) - E{f) = Ed + U>0 ^ Ea + U/2> -[7/2, 

or (Fig. do]). 

U/2 > Ed + U/2 > -U/2. 



(23) 



(24) 



Under these conditions, a local moment is well-defined provided the temperature is smaller 
than the valence fluctuation scale Typ = max(£'j + U,—Ef). At lower temperatures, the 
atom behaves exclusively quantum top. 

Virtual bound-state formation. When the magnetic ion is immersed in a sea of electrons, 
the f-electroi is within the core of the a tom hybridize with the Bloch states of surrounding 



electron sea (jBlandin and Friedell . 



(25) 



19581 ) to form a resonance described by 

Hresonance = ^ ^kUka + ^ \^V{k)cl^f^ + V (k)* f^Cka 

where the hybridization matrix element V{k) = (k| V^tomid/) is the overlap of the atomic 
potential between a localized f-state and a Bloch wave. In the absence of any interactions, 
the hybridization broadens the localized f-state, producing a resonance of width 

A = 7r^|y(A:)|25(ek-/i) =vryV (26) 
k 

where V'^ is the average of the hybridization around the Fermi surface. 

There are two complimentary ways to approach the physics of the Anderson model 

- the "atomic picture", which starts with the interacting, but isolated atom {V{k) = 0), and 
considers the effect of immersing it in an electron sea by slowly dialing up the hybridization. 
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- the "adiabatic picture" which starts with the non- interacting resonant ground-state {U = 0), 
and then considers the effect of diahng up the interaction term U. 

These approaches paint a contrasting and, at first sight, contradictory picture of a local moment in 
a Fermi sea. From the adiabatic perspective, the ground-state is always a Fermi liquid (see II.B.2|) . 
but from atomic perspective, provided the hybridization is smaller than U one expects a local 
magnetic moment, whose low lying degrees of freedom are purely rotational. How do we resolve 
this paradox? 

Anderson's original work provided a mean-field treatment of the interaction. He found that 
at interactions larger than Uc ~ vrA local moments develop with a finite magnetization M = 
(n|) — (nj). The mean field theory provides an approximate guide to the conditions required 
for moment formation, but it doesn't account for the restoration of the singlet symmetry of the 
ground-state at low temperatures. The resolution of the adiabatic and the atomic picture derives 
from quantum spin fluctuations which cause the local moment to "tunnel" on a slow time-scale Tsf 
between the two degenerate "up" and "down" configurations. 



ej + fl^e- + fl 



(27) 



These fluctuations are the origin of the Kondo effect. From the energy uncertainty principle, 
below a temperature Tk at which the thermal excitation energy ksT is of order the characteristic 
tunneling rate — , a paramagnetic state with a Fermi liquid resonance will form. The characteristic 



width of the resonance is then determined by the Kondo energy knTx ^ — 



h rp 



this resonance was first deduced by Abrikosov and Suhl (|Abrikosovl . 



1965 



Suhi 



he existence of 



19651 ). but it is 



more frequently called the "Kondo resonance" . From perturbative renormaliz ation group reasoning 



1983 


Wieemann . 


1980a 


b) 



19 781) a nd the Bethe ansatz solution of the Anderson model ( Okiii and Kawakamil . 

A, the Kondo scale depends exponentially 
-C//2, 



on U . In the symmetric Anderson model, where Ej 

Tk = 



2UA 
— — exp 



ttU 
'8A 



(28) 



The temperature Tk marks the crossover from a a high temperature Curie law x ^ t susceptibility 
to a low temperature paramagnetic susceptibility x ~ ^/Tk- 



23 



2. Adiabaticity and the Kondo resonance 



A central quantity in the physics of f-electron systems is the f-spectral function, 

Af{uj) = -ImGfiuj - i6) 



(29) 



where Gf{uj) = —iJ^^dt{Tfcr{t)fa{0))e^'^^ is the Fourier transform of the time-ordered f-Green's 
function. When an f-electron is added, or removed from the f-state, the final state has a distribution 
of energies described by the f-spectral function. From a spectral decomposition of the f-Green's 
function, the positive energy part of the f-spectral function determines the energy distribution for 
electron addition, while the negative energy part measures the energy distribution of for electron 



removal: 



' Energy distribution of state formed by adding one f-clectron. 



Afiu;) = < 



5^|(A|/t|0o>| 6{u;-[Ex-Eo]), {oo > 0) 

A 

Y,\iMfM\^ S{u; - [Eo - Ex]), (u; < 0) 



(30) 



Energy distribution of state formed by removing an f-elcctron 



where £"0 is the energy of the ground-state, and Ex is the energy of an excited state A, formed 
by adding or removing an f-electron. For negative energies, this spectrum can be measured by 
measuring the energy distribution of photo-electrons produced by X-ray photo-emission, while 
for positive energies, th e spectral function can be measured from inverse X-ray photo-emission 



(Allen et al. 



1986 



19831 ). The weight beneath the Fermi energy peak, determines the f-charge of 



the ion 



(31) 



In a magnetic ion, such as a Cerium atom in a 4/^ state, this quantity is just a little below unity. 

Fig. ()lll ) illustrates the effect of the interaction on the f-spectral function. In the non- 
interacting limit {U = 0), the f-spectral function is a Lorentzian of width A. If we turn on the 
interaction U , being careful to shifting the f-level position beneath the Fermi energy to maintain 
a constant occupancy, the resonance splits into three peaks, two at energies uj = Ef and lo = 
Ef + U corresponding to the energies for a valence fluctuation, plus an additional central "Kondo 
resonance" associated with the spin-fluctuations of the local moment. 

At first sight, once the interaction is much larger than the hybridization width A, one might 
expect there to be no spectral weight left at low energies. But this violates the idea of adiabaticity. 
In fact, there are always certain adiabatic invariants that do not change, despite the interaction. 
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Af((0) 




0) = E^+U 



FIG. 11 Schematic illustrating the evaluation of the f-spectral function Af{u!) as interaction strength U is 
turned on continuously, maintaining a constant f-occupancy by shifting the bare f-level position beneath 
the Fermi energy. The lower part of diagram is the density plot of f-spectral function, showing how the 
non-interacting resonance at U — splits into an upper and lower atomic peak at ui — Ef and uj — Ef + U. 



One such quantity is the phase shift 5f associated with the scattering of conduction electrons off 
the ion; another is the height of the f-spectral function at zero energy, and it tu rns out that these 
two quantities are related. A rigorous result due to Langreth (jLangrethl . 



19661 ). tells us that the 



spectral function at w = is diretly determined by the f-phase shift, so that its non-interacting 
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value 



Af{Lo = 0) 



sin^ 5 j 
vrA 



(32) 



is preserved by adiabaticity. Langreth's result can be heuristically derived by noting that (5/ is the 
phase of the f-Green's function at the Fermi energy, so that G f{Q — i5)~^ = ||Gj^(0)|e~*'^^. Now in 
a Fermi liquid, the scattering at the Fermi energy is purely elastic, and this implies that ImGj^{0 — 
iS) = A, the bare hybridization width. From this it follows that /mGj^(O) = |Gj^(0)| sinJj = A, 
so that Gf{0) = e^^f /{Asm6), and the above result follows. 

The phase shift Sf is set via the Friedel sum rule, according to which the sum of the up and 
down scattering phase shifts, gives the total number of f-bound-electrons, or 



Uf. 



(33) 



for a two-fold degenerate /— state. At large distances, the wavefunction of scattered electrons 
ipf{r) ~ sm{kpr + 6f)/r is "shifted inwards" by by a distance Si/kp = (Air/2) x (di/ir). This 
sum rule is sometimes called a "node counting" rule, because if you think about a large sphere 
enclosing the impurity, then each time the phase shift passes through vr, a node crosses the spherical 
boundary and one more electron per channel is bound beneath the Fermi sea. Friedel's sum rule 



holds for interactin g electrons, providing the groun c 



interacting system (ILanger and Ambegaokar 



1961 



state is adiabat ically accessible from the non- 



Langreth 



19661 ). Since = 1 in an state. 



the Friedel sum rule tells us that the phase shift is 7r/2 for a two-fold degenerate /— state. In other 
words, adiabaticity tell us that the electron is resonantly scattered by the quenched local moment. 
Photo-emission studies do reveal the three-peaked st ructure cha r acteri stic of the Anderson 



19831 ) (see Fig. O. Ma- 



model in many Ce systems, such as Celr2 and CeRu2 (jAllen et al. 
terials in which the Kondo resonance is wide enough to be resolved are more "mixed valent" 
materials in which the f- valence departs significantly from unity. T hree peaked struc tures have 



also been observed in certain U 5f materials such as U Pt^ and U AI2 (jAllen et al 



19851 ) materials , 



but it has not yet been resolved in UBei^. A three peaked structure has recently been observed 
in 4f Yb materials, such as YbPd^, where the 4/^^ configur ation contains a single / hole, so that 



the positions of the three peaks are reversed relative to Ce (jLiu et al. 



1992 ). 
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FIG. 12 Showing spectral functions for three different Cerium f-electron materials, measured using X-ray 
photoemission (be low the Fermi energy ) and inverse X-ray photoemission (above the Fermi energy) after 



(Allen et al. 



1983|). CeAl is an antiferromagnet and does not display a Kondo resonance. 



B. Hierachies of energy scales 



1. Renormalization Concept 



To understand how a Fermi liquid emerges when a local moment is immersed in a quantum sea of 
electrons, theorists had to connect physics on on several widely spaced energy scales. Photoemission 
shows that the characteristic energy to produce a valence fluctuation is of the order of volts, or 
tens of thousands of Kelvin, yet the the characteristic physics we are interested in occurs at scales 
hundreds, or thousands of times smaller. How can we distill the essential effects of the atomic 
physics at electron volt scales on the low energy physics at millivolt scales? 



The essential tool for this task is the " r enormalization group" (lAnderson 



Anderson and Yuval 



196S 



1970!, 



1971 



Nozieres 



1976 



Nozieres and Blandin 



1980 



197C 



Wilson 



19731 : 



1973), 



based on the idea that the physics at low energy scales only depends on a small subset of "relevant" 
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FIG. 13 (a) Cross-over energy scales for the Anderson model. At scales below A/, valence fluctuations into 
the doubly occupied state are suppressed. All lower energy physics is described by the infinite U Anderson 
model. Below A/j, all valence fluctuations are suppressed, and the physics involves purely the spin degrees of 
freedom of the ion, coupled to the conduction sea via the Kondo interaction. The Kondo scale renormalizes 
to strong coupling below A//j , and the local moment becomes screened to form a local Fermi liquid, (b) 
Illustrating the idea of renormalization group flows towards a Fermi liquid flxed point. 

variables from the original microscopic Hamiltonian. The extraction of these relevant variables is 
accomplished by "renormalizing" the Hamiltonian by systematically eliminating the high energy 
virtual excitations and adjusting the low energy Hamiltonian to take care of the interactions that 
these virtual excitations induce in the low energy Hilbert space. This leads to a family of Hamilto- 
nian's H{A), each with a different high-energy cut-off A, which share the same low energy physics. 

The systematic passage from a Hamiltonian H{A) to a renormalized Hamiltonian H{A') with 
a smaller cutoff A' = A/6 is accomplished by dividing the the eigenstates of H into a a low energy 
subspace {L} and a high energy subspace {H}, with energies |e| < A' = A/b and a |e| G [A', A] 
respectively. The Hamiltonian is then broken up into terms that are block-diagonal in these 
subspaces, 



H 





yt 1 


_ V 





(34) 



where V and V'^ provide the matrix elements between {L} and {-ff}. The effects of the V are then 
taken into account by carrying out a unitary (canonical) transformation that block-diagonalizes 
the Hamiltonian, 



H{A) UH{A)U^ 











Hh. 



(35) 
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The renormalized Hamiltonian is then given by H{A') = Hl = Hl+SH . The flow of key parameters 
in the Hamiltonian resulting from this process is called a renormalization group flow. 

At certain important cross-over energy scales, large tracts of the Hilbert space associated with 
the Hamiltonian are projected out by the renormalization process, and the character of the Hamil- 
tonian changes qualitatively. In the Anderson model, there are three such important energy 
scales. (fT3]) 



• Kj = Ef + U where valence fluctuations e + ^ f"^ into the doubly occupied state are 
eliminated. For A << A/, the physics is described by the infinite U Anderson model 

V{k) cl^Xoa + X^ocka l+Ef^X^^, (36) 

k,cr k,a o" 

where X,, = : (j){f^ : a\, Xq. = and X^o = \ f ■ are "Hubbard 

operators" that connect the states in the projected Hilbert space with no double occupancy. 

• A// ~ \Ef\ = —Ef, where valence fiuctuations into the empty state ^ Z*^ + e~ are 
eliminated to form a local moment. Physics below this scale is described by the Kondo 
model 

• A = Tk, the Kondo temperature below which the local moment is screened to form a 
resonantly scattering local Fermi liquid. 

In the symmetric Anderson model, A/ = A//, and the transition to local moment behavior occurs 
in a one-step crossover process. 



2. Schrieffer Wolff transformation 



The unitary, or canonical transformation tha t eliminates the charge fluctuations at scales A j and 



1969 



Schrieffer and Wolfil . 



A. J J w as first carried out by Schrieffer and Wolff (jCoqblin and Schrieffe: 
19661 ). who showed how this model gives rise to a residual antiferromagnetic interaction between 
the local moment and conduction electrons. The emergence of this antiferromagnetic interaction is 
associated with a process called "superexchange" : the virtual process in which an electron or hole 
briefly migrates off the ion, to be immediately replaced by another with a different spin. When 
these processes are removed by the canonical transformation, they induce an antiferromagnetic 
interaction between the local moment and the conduction electrons. This can be seen by considering 
the two possible spin exchange processes 

ef + fl ^ el + /t A^, ~ [/ + E^ 
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'^11 



-E 



f 



(37) 



Both process requires that the f-electron and incoming particle are in a spin-singlet. From second 
order perturbation theory, the energy of the singlet is lowered by an amount —2 J where 



J = V' 



1 1 

+ 



(38) 



AEi AE2 

and the factor of two derives from the two ways a singlet can emit an electron or hole into the 
continuum ^ and V ~ V{kp) is the hybridization matrix element near the Fermi surface. For the 
symmetric Anderson model, where AEi = AEu = U /2, J = ^V"^ /U . 

If we introduce the electron spin density operator (t(0) = ^, c^Q^a/jCfc'/j, where J\f is the 

number of sites in the lattice, then the effective interaction will have the form 



Hk = -2JPs= 



(39) 



where Ps= 



1 1 



(t(0) • Sf is the singlet projection operator. If we drop the constant term, 



then the effective interaction induced by the virtual charge fluctuations must have the form 

Hk = Ja{0) • Sf 



3. The Kondo Effect 



(40) 



where Sf is the spin of the localized moment. The complete "Kondo Model", H = He + Hk 
describing the conduction electrons and their interaction with the local moment is 



(41) 



The anti-ferromagnetic sign of the super-exchange interaction J in the Kondo Hamiltonian is 
the origin of the spin screening physics of the Kondo effect. The bare interaction is weak, but 
the spin fluctuations it induces have the effect of antiscreening the interaction at low energies, 



renormalizing it to larger and larger va 



scaling procedure ([Anderson 



1971 



ues. To see this, we follow a Anderson's "Poor Man's" 



I973I ). which takes advantage of the observation that at small 



To calculate the matrix elements associated with valence fluctuations, take 



and 



cljO) 



then (c^l VctUlfc^) = V2V and {f\ VfU.lfc^) = V2V 
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J the renormalization in the Hamiltonian associated with the block-diagonahzation process 5H = 
Hl — Hl is given by second-order perturbation theory: 



SH^ = {a\6H\b) = i [T,b{Ea) + TabiE^)] 



(42) 



where 



\A)e{H} 



lo-Ea 



(43) 



is the many body "t-matrix" associated with virtual transitions into the high-energy subspace {H}. 
For the Kondo model, 



V = VhJS{Q) ■ SoPl 



(44) 



where Vh projects the intermediate state into the high energy subspace while Vl projects the initial 
state into the low energy subspace. There are two virtual scattering processes that contribute the 
the antiscreening effect, involving a high energy electron (I) or a high energy hole (II). 
Process I is denoted by the diagram 




and starts in state \b) = \ka, a) , passes through a virtual state |A) = \cl,„aa") where eu" lies at high 
energies in the range e^" G [A/6, A] and ends in state \a) = \k'P,a'). The resulting renormalization 

1 



{k'(3,a'\T'iE)\ka,a)= = ^ 

ej.//e[A-(5A,A] 



J^pSA 



1 



E-A 



E — ek" 



(45) 



In Process II , denoted by 




the formation of a virtual hole excitation |A) = Ck"x\a-") introduces an electron line that crosses 
itself, introducing a negative sign into the scattering amplitude. The spin operators of the conduc- 
tion sea and antiferromagnet reverse their relative order in process II, which introduces a relative 
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minus sign into the T-matrix for scattering into a high-energy hole-state, 

1 



{k'Pa'\T^^^\E)\k 



aa) 



E 

efc//e[-A,-A+<5A] 



1 



E - {ek + ey - Cfc") 



a Qb\ 



E-A 

where we have assumed that the energies and e^/ are neghgible compared with A. 
Adding (Eq. HS]) and (Eq. HI]) gives 

,1 , ^jj _ -PpSA,, 



(46) 



s: Tjint 



T +T 
A 



A 



(47) 



so the high energy virtual spin fluctuations enhance or "anti-screen" the Kondo coupling constant 



J(A') = J(A) + 2J> 



2 SA 



A 



(48) 



If we introduce the coupling constant g = pJ, recognizing that din A = — ^we see that it satisfies 



dg 
dlnA 



(3ig)=-2g^ + Oig'). (49) 



This is an example of a negative (3 function: a signature of an interaction which grows with 
the renormalization process. At high energies, the weakly coupled local moment is said to be 
asymptotically free. The solution to the scaling equation is 



9(A') 



9o 



l-2goHA/A') 



(50) 



and if we introduce the "Kondo temperature" 

Tk = D exp 

we see that this can be written 

25(A'; 



1 ■ 



(51) 



1 



(52) 



ln(A/r,^) 

so that once A' ~ Tr-, the coupling constant becomes of order one - at lower energies, one reaches 
"strong coupling" where the Kondo coupling can no longer be treated as a weak perturbation. 
One of the fascinating things about this flow to strong coupling, is that in the limit Tr- « D, all 
explicit dependence on the bandwidth D disappears and the Kondo temperature Tk is the only 
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intrinsic energy scale in the physics. Any physical quantity must depend in a universal way on 
ratios of energy to Tr-, thus the universal part of the Free energy must have the form 



F(r) = TkHt/Tk] 



(53) 



where is universal. We can also understand the resistance created by spin- flip scattering 



off a magnetic impurity in the same way. The resistivity is given by pi 
scattering rate must also have a scaling form 



r(T,H) = -^^2{T/Tk,H/Tk) 
P 



-t(T,H) where the 



(54) 



where p is the density of states (per spin) of electrons and rii is the concentration of magnetic 
impurities and the function <I>2(t,/i) is universal. To leading order in the Born approximation, 
the scattering rate is given by r = 2-Kpj'^S{S + 1). = ^^'^(^+^) (gg)^ where = ^(Ao) is the bare 
coupling at the energy scale that moments form. We can obtain the behavior at a finite temperature 
by replacing qq — > g{K = 2itT), where upon 

2ttS{S + 1) 1 



t{T) 



(55) 



P A\xi^{2ttT/Tk) 

gives the leading high temperature growth of the resistance associated with the Kondo effect. 

The kind of perturbative analysis we have gone through here takes us down to the Kondo 
temperature. The physics at lower energies requires corresponds to the strong coupling limit of the 
Kondo model. Qualitatively, once Jp » 1, the local moment is bound into a spin-singlet with a 
conduction electron. The number of bound-electrons is = 1, so that by the Friedel sum rule 
(eq. [33]) in a paramagnet the phase shift d-^ = = vr/2, the unitary limit of scattering. For more 
details about the local Fermi liquid th at forms, we refer the reader to the accompanying chapter 



on the Kondo effect by Barbara Jones ((Jones 



2003) 



4-. Doniach's Kondo Lattice Concept 



The discovery of heavy electron metals prompted Doniach (|Doniachl . 119771 ) to make the radical 
proposal that heavy electron materials de rive from a d ense lattice version of the Kondo effect. 



described by the Kondo Lattice model (jKasuyal . 



19561) 



(56) 



ko- j 

In effect, Doniach was implicitly proposing that the key physics of heavy electron materials resides 
in the interaction of neutral local moments with a charged conduction electron sea. 
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FIG. 14 Spin polarization around magnetic impurity contains Friedel oscillations and induces an RKKY 
interaction between the spins 

Most local moment systems develop antiferromagnetic order at low temperatures. A magnetic 
moment at location xq induces a wave of "Friedel" oscillations in the electron spin density (Fig. 



((7(x)) = - Jx(x - xo)(S'(xo)) 



where 



k,fc' 



/(gk) - /(gk') \ gi(k-k')-x 

ek' - ek 



(57) 



(58) 



is the non-local susceptibility of the metal. The sharp discontinuity in the occupancies /(e^) at 
the Fermi surface is responsible for Friedel oscillations in induced spin density that decay with a 
power-law 

cos 2kFr 



{a{r)) Jp- 



(59) 



where p is the conduction electron density of states and r is the distance from the impurity. If a 
second local moment is introduced at loc ation x, it couples to this Friedel oscillation with energy 



J(5(x)-a(x)) giving rise to the "RKKY" (IKasuva 
magnetic interaction, 



1953; 



Ruderman and Kittel 



1950; 



Yosidal . 



19571 ^ 



Jrkky{x-x') 



where 



Hrkky = -J'x(x - x') S^(x) • 5(x'). 



2 cos2kFr 

jRKKY{r) ~ -J P . 



(60) 



(61) 



In alloys containing a dilute concentration of magnetic transition metal ions, the oscillatory RKKY 
interaction gives rise to a frustrated, glassy magnetic state known as a "spin glass". In dense 
systems, the RKKY interaction typically gives rise to an ordered antiferromagnetic state with a 
Neel temperature T/v of order J^p. Heavy electron metals narrowly escape this fate. 
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Tk Dp.Tp[—l/Jp] 



Jp 




RKKY 



T^>T 



RKKY 



FIG. 15 Doniach diagram, illustrating the antiferromagnetic regime, where Tr- < T^kky and the heavy 
fermion regime, where Tk > Trkky- Experiment has told us in recent times that the transition between 
these two regimes is a quantum critical point. The effective Fermi temperature of the heavy Fermi liquid is 
indicated as a solid line. Circumstantial experimental evidence suggests that this scale drops to zero at the 
antiferromagnetic quantum critical point, but this is still a matter of controversy. 

Doniach argued that there are two scales in the Kondo lattice, the single ion Kondo temperature 
Tk and Trkky, given by 



When Jp is small, then TjixKY is the largest scale and an antiferromagnetic state is formed, but 
when the Jp is large, the Kondo temperature is the largest scale so a dense Kondo lattice ground- 
state becomes stable. In this paramagnetic state, each site resonantly scatters electrons with a 
phase shift ~ tt/2. Bloch's theorem then insures that the resonant elastic scattering at each site 
will act coherently, forming a renormalized band of width ~ Tk (Fig. [T5]) . 

As in the impurity model, one can identify the Kondo lattice ground-state with the large U 
limit of the Anderson lattice model. By appealing to adiabaticity, one can then link the excitations 
to the small U Anderson lattice model. According to this line of argument, the quasiparticle Fermi 



Tk = Z)e-i/(2Jp) 
Trkky = J^P 



(62) 
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Screening Cloud 




Overlap? 



Exhaustion 
N(0)Tk « 1 ? 




Composite Heavy 
Femiion 



New states injected 
into Fermi sea. 



FIG. 16 Contrasting (a) tlie "screening cloud" picture of tlie Kondo effect with (b) the composite fermion 
picture. In (a), low energy electrons form the Kondo singlet, leading to the exhaustion problem. In (b) the 
composite heavy electron is a highly localized bound-state between local moments and high energy electrons 
which injects new electronic states into the conduction sea at the chemical potential. Hybridization of these 
states with conduction electrons produces a singlet ground-state, forming a Kondo resonance in the single 
impurity, and a coherent heavy electron band in the Kondo lattice. 



surface volume must count the number of conduction and f-electrons (iMartinl . Il982l ) even in the 
large U limit, where it corresponds to the number of electrons plus the number of spins 



(27r)3 



rie + n. 



spms' 



(63) 



Using top ology, and certa in basic assumptions about the response of a Fermi liquid to a flux, 



Oshikawa (Oshikaw: 



2000) has been able to short-circuit this tortuous path of reasoning, proving 



that the Luttinger relationship holds for the Kondo lattice model vi^ithout reference to its finite U 
origins. 

There are however, aspects to the Doniach argument that leave cause for concern: 

• it is purely a comparison of energy scales and does not provide a detailed mechanism con- 
necting the heavy fermion phase to the local moment antiferromagnet. 
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simple estimates of the value of Jp required for heavy electron behavior give an artificially 
large value of the coupling constant Jp ~ 1. This issue was later resolved by the observation 
that large spin degeneracy 2j + l of the spin-orbit coupled moments, which can be as large as 
= 8 in Yh materials , enha nces the rate of scaling to strong coupling, leading to a Kondo 



temperature (jColemanl . 



1983) 



Tk = D{NJp)N exp 



NJp 



(64) 



Since the scaling enhancem ent effect stretches out across decades of energy, it is largely 



robust against crystal fields (jMekata et al 



19861). 



Nozieres' exhaustion paradox (jNozieres 



19851 ). If one considers each local moment to be 



magnetically screened by a cloud of low energy electrons within an energy Tk of the Fermi 
energy one arrives at an "exhaustion paraodox". In this interpretation, the number of 
electrons available to screen each local moment is of order Tk/D << 1 per unit cell. Once 
the concentration of magnetic impurities exceeds ^ ~ 0.1% for {Tk = l^K, D = 10^-fC), 
the supply of screening electrons would be exhausted, logically excluding any sort of dense 
Kondo effect. Experimentally, features of single ion Kondo behavior persist to much higher 
densities. The resolution to the exhaustion paradox lies in the more modern perception that 
spin-screening of local moments extends up in energy , from the Kondo scale Tk out to the 
bandwidth. In this respect, Kondo screening is reminiscent of Cooper pair formation, which 
involves electron states that extend upwards from the gap energy to the Debye cut-off. From 
this perspective, t he Kondo length sca le ^ ~ vf/Tk is analogous to the coherence length of 



a superconductor (|Burdin et al. 



20001 ). defining the length scale over which the conduction 



spin and local moment magnetization are coherent without setting any limit on the degree 
to which the correlation clouds can overlap ( Fig. [T6|) . 

C. The Large N Kondo Lattice 

1. Gauge theories, Large N and strong correlation. 



The "standard model" for metals is built upon the expansion to high orders in the strength 
of the interaction. This approach, pioneered by Landau, and later formulated in the language 



of finite temperature perturbation 



(jLandaul . 



19571 : 



theory 



Luttinger and Ward 



19601 : 



jy Pitaevksii, Luttinge r ^ Ward, Nozieres and others 



Nozieres and Luttinged . 



1962 



Pitaevskii 



196Q!), pro- 



vides the foundation for our understanding of metallic behavior in most conventional metals. 
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The development of a parallel formalism and approach for strongly correlated electron systems 
is still in its infancy, and there is no universally accepted approach. At the heart of the problem 
are the large interactions which effectively remove large tracts of Hilbert space and impose strong 
constraints on the low-energy electronic dynamics. One way to describe these highly constrained 
Hilbert spaces, is through the use of gauge theories. When written as a field theory, local con- 
straints manifest themselves as locally conserved quantities. General principles link these conserved 
quantities with a set of gauge symmetries. For example, in the Kondo lattice, if a spin S = 1/2 
operator is represented by fermions, 

S3=fli^^ /.73, (65) 



'i - Jja j ^•'J'^' 

then the representation must be supplemented by the constraint nf{j) = 1 on the conserved f- 
number at each site. This constraint means one can change the phase of each f-fermion at each 
site arbitrarily 

fj^e^'^^fj, (66) 
without changing the spin operator Sj or the Hamiltonian. This is the local gauge symmetry. 



Similar issues also arise in the infinite U Anderson or Hubbard models where t 



re "no double 



occupancy" co nstraint can be established by using a slave boson representation (jBarnesl . 



Colemanl . 



19761 : 



19841 ) of Hubbard operators: 

XMj) = flbj, Xo^ij) = b]f,, (67) 

where /j^ creates a singly occupied f-state, = 1/"*^;^) while 6^ creates an empty state, 

6^10) = In the slave boson, the gauge charges 

Q. = E4/i- + % (68) 

cr 

are conserved and the physical Hilbert space corresponds to Qj = 1 at each site. The gauge 
symmetry is now /jo- — > e^^^ fja, bj e^^^hj. These two examples illustrate the link between strong 
correlation and gauge theories. 

strong correlation ^ constrained Hilbert Space <-> gauge theories (69) 

A key feature of these gauge theories, is the appearance of "fractionalized fields" which carry either 
spin or charge, but not both. How then, can a Landau Fermi liquid emerge within a Gauge theory 
with fractional excitations ? 
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Some have suggested that Fermi hquids can not reconstitute themselves in such strongly con- 
strained gauge theories. Others have advocated against gauge theories, arguing that the only 
reliable way forward is to return to "real world" models with a full fermionic Hilbert space and 
a finite interaction strength. A third possibility is that the gauge theory approach is valid, but 
that heavy quasiparticles emerge as bound-states of gauge particles. Quite independently of one's 
position on the importance of gauge theory approaches, the Kondo lattice poses a severe compu- 
tational challenge, in no small part, because of the absence of any small parameter for resumed 
perturbation theory. Perturbation theory in the Kondo coupling constant J always fails below the 
Kondo temperature. How then, can one develop a controlled computational tool to explore the 
transition from local moment magnetism to the heavy fermi liquid? 

One route forward is to seek a family of models that interpolates between the models of physical 
interest, and a limit where the physics can be solved exactly. One approach, as we shall discuss 
later, is to consider Kondo lattices in variable dimensions d, and exp and in powers of 1/d about 



the limit of infinite dimensionality (jGeorges et al. 



1996 



Jarrelll . 



19951 ). In this limit, electron self- 



energies become momentum independent, the basis of the dynamical mean- field theory. Another 
approach, with the advantage that it can be married with gauge theory, is the use of large 
expansions. The idea here is to generalize the problem to a family of models in which the f-spin 
degeneracy A = 2j + 1 is artificially driven to infinity. In this extreme limit, the key physics is 
captured as a mean-field theory, and finite A properties are obtained through an expansion in the 
small parameter 1/A. Such large A^ expansi ons have played an irn portant role in the co ntext of 



the s pherical model of statistical mechanics (iBerlin and Kad . 



1952) and in field theory (jWittenl . 



19781 ). The next section discusses how the gauge theory of the Kondo lattice model can be treated 



in a large A expansion. 



2. Mean field theory of the Kondo lattice 



Quantum large A expansions are a kind of semi-classical limit where 1 /N ~ h plays the role of 
a synthetic Planck's constant. In a Feynman path integral 



{xf{t)\xi,0) = / Vlx] exp 



(70) 



where S is the classical action and the quantum action A = is "extensive" in the variable ^. 
When i — > oo, fluctuations around the classical trajectory vanish and the transition amplitude is 



entirely determined by the classical action to go from i to /. A large A expansion for the partition 
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function Z of a quantum system involves a path integral in imaginary time over the fields <f) 

Z = [ V[4>]e-^^^'f''^^ (71) 



where NS is the action (or free energy) associated with the field configuration in space and time. 
By comparison, we see that the large N limit of quantum systems corresponds to an alternative 
classical mechanics where 1/N ^ fi emulates Planck's constant and new types of collective behavior 
not pertinent to strongly interacting electron systems, start to appear. 
Our model for a Kondo lattice of spins localized at sites j is 

ko- j 

where 

Hi{j) = j^Sa(3{3)c]pCjc (73) 



19691 ) 



is the Coqblin Schrieffer form of the Kondo interaction Hamiltonian (iCoqblin and Schrieffer , 
between an f-spin with N = 2j + \ spin components and the conduction sea. The spin of the local 
moment at site j is represented as a bilinear of Abrikosov pseudo-fermions 

So.p{j) = flhp - (74) 

and 

=^Vc| e-'^-^^ (75) 

creates an electron localized at site j, where Ns is the number of sites. 

Although this is a theorists' idealization - a "spherical cow approximation", it nevertheless 
captures key aspects of the physics. This model ascribes a spin degeneracy of = 2j + 1 to the 
both the f-electrons and the conduction electrons. While this is justified for a single impurity, a 
more realistic lattice model requires the introduction of Clebsch Gordon coefficients to link the 
spin-1/2 conduction electrons with the spin-j conduction electrons. 

To obtain a mean- field theory, each term in the Hamiltonian must scale as N . Since the 
interaction contains two sums over the spin variables, this criterion is met by rescaling the coupling 
constant replacing J — > Another important aspect to this model, is the constraint on charge 
fluctuations, which in the Kondo limit imposes the constraint nj = 1. Such a constraint can be 
imposed in a path integral with a Langrange multiplier term A(n/ — 1). However, with nj = 1, 
this is not extensive in N, and can not be treated using a mean-field value for A. The resolution is 
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to generalize the constraint to nj = Q, where Q is an integer chosen so that as grows, q = Q/N 
remains fixed. Thus, for instance, if we are interested in = 2, this corresponds to q = nj/N = i. 
In the large limit, it is then sufficient to apply the constraint on the average (uf) = Q, through 
a static Lagrange multiplier coupled to the difference (nj — Q). 

The next step is to carry out a "Hubbard Stratonovich" transformation on the interaction 



J 



HiU) = -j^{cU,p){flc 



(76) 



Here we have absorbed the term —j^njc^-Cja derived from the spin-diagonal part of (j74p by a shift 



in the chemical potential. This interaction has the form —gA^A, with g 



N 



and 



A = fj^Cja, which we factorize using a Hubbard Stratonovich transformation. 



jA^A A^V + VA + 



vv 



(77) 



so that (jLacroix and Cvrot 



1981 



Read and Newns 



1983a|) 



Hi{j) ^ Hi[V,j] = Vj {c]j,J + /t c, J Vj + N 



J 



(78) 



This is an exact transformation, provided the Vj{T) are treated as ffuctuating variables inside a 
path integral. The Vj can be regarded as a spinl ess exchange boso n for the Kondo effect. In the 



parallel treatment of the infinite Anderson model (|Colemanl . 
field associated with valence ffuctuations. 
In diagrams: 



1987al ). Vj = Vbj is the "slave boson" 
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where 



(81) 



kcr 



1986 



Read and Newns 



This i s the "Read Newns" path integral formulation (jAuerbach and Levin 
1983al ) of the Kondo lattice model. The path integral contains an outer integral J !)[]/, \] over 
the gauge fields Vj and Aj(r), and an inner integral J T>[c, f] over the fermion fields moving in the 
environment of the gauge fields. The inner path integral is equal to a trace over the time-ordered 
exponential of H[V, A]. 

Since the action in this path integral grows extensively with N, the large N limit is saturated 
by the saddle point configurations of V and A, eliminating the the outer integral in (j8ip . We 
seek a translationally invariant, static, saddle point where Aj(r) = A and V^(t) = V. Since the 
Hamiltonian is static, the interior path integral can be written as the trace over the Hamiltonian 
evaluated at the saddle point, 



Tre 



-PHmft 



{N oo) 



(82) 



where 




The saddle point is determined by the condition that the Free energy F = —T In Z is stationary 
with respect to variations in V and A. To impose this condition we need to diagonalize Hmft and 
compute the Free energy. First we rewrite the mean field Hamiltonian in momentum space. 



ek 


V 


/ Ck<7 \ 


V 


A 





Xq 



(84) 



where 



/, 



fccr 



fja 



(85) 



is the Fourier transform of the /—electron field. This Hamiltonian can then be diagonalized in the 
form 



Hmft = (aL,&^ 



/ ^ \^"ko-' "kcr 
kfT 



^k+ 

Ek- 




J 



Xq 



(86) 
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FIG. 17 (a) Dispersion produced by the injection of a composite fermion into the conduction sea. (b) 
Renormahzed density of states, showing "hybridization gap" (A^). 

where a\^^ and b^^^ are linear combinations of c|^^ and ft which describe the quasiparticles of the 
theory. The momentum state eigenvalues E = are the roots of the equation 

{E-e^,){E-\)-\V\'' = 



Det 















1 











(87) 



so 



Ek - A 



(88) 



are the energies of the upper and lower bands. The dispersion described by these energies is shown 
in Fig. dZ] . Notice that: 

• hybridization between the f-electron states and the conduction electrons builds an upper 
and lower Fermi band, separated by an indirect "hybridization gap" of width = Eg[+) — 
Eg{—) ~ Tk, where 

y2 



Eg{±)=\± 



(89) 



and ^D± are the top and bottom of the conduction band. The "direct" gap between the 
upper and lower bands is 2\V\. 

From (j87p . the relationship between the energy of the heavy electrons {E) and the energy of 
the conduction electrons (e) is given by e = E — \ V'\^ /{E — X), so that the density of heavy 
electron states p*{E) = '^\^±^{E — E^'^) is related to the conduction electron density of 
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states p(e) by 



de . ^ / IV'P \ I p ( 1 + (J^ A)^ ) outside hybridization gap, 



p*(i?) = p— = p(e) 1 + 



tfE ' V - A)2 



inside hybridization gap, 

(90) 



so the "hybridization gap" is flanked by two sharp peaks of approximate width Tk- 

• The Fermi surface volume expands in response to the injection of heavy electrons into the 
conduction sea, 

^«''7^ = (7^E^k.) = Q + ne (91) 

where is the unit cell volume, nko- = ako-^kcr + ^ko-^ko- is the quasiparticle number operator 
and ric is the number of conduction electrons per unit cell. More instructively, if rif, = ndoP 
is the electron density, 

e~ density q.particle density 

^ _ 9 

(27r)3 a 

+ve background 



so the electron density ric divides into a contribution carried by the enlarged Fermi sea, whose 
enlargement is compensated by the development of a positively charged background. Loosely 
speaking, each neutral spin in the Kondo lattice has "ionized" to produce Q negatively 
charged heavy fermions, leaving behind a Kondo singlet of charge +Qe (Fig. [TSl). 



To obtain V and A we must compute the Free energy 

1 + e-^^i-i 



- = -Tyin 



+ (^-Ag). (93) 



At T = 0, the Free energy converges the ground-state energy E'oi given by 



En 



Using (f90l ) , the total energy is 



(94) 



E /"O p /ii/p 



D^p A /Ae , , 

-^+-lnf^)-A, (95) 
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Heavy Electrons 



Kondo singlets: 
charged backgrouric 

FIG. 18 Schematic diagram from (|Coleman et al\ . \2005a ).(a) High temperature state: smaU Fermi surface 



with a background of spins; (b)Low temperature state where large Fermi surface develops against a back- 
ground of positive charge. Each spin "ionizes" into Q heavy electrons, leaving behind a a Kondo singlet 
with charge +Qe 

where we have assumed that the upper band is empty and the lower band is partially filled. 

i_ 

Tk = De as before. The first term in (j95l ) is the conduction electron contribution to the 
energy Ec/Nug, while the second term is the lattice "Kondo" energy EK/Nj\f^. If now we impose 
the constraint = (nj) — Q = then A = ^ so that the ground-state energy can be written 



Ek ^ a / Ae \ 
This energy functional has a "Mexican Hat" form, with a minimum at 



(96) 



(97) 



confirming that A ~ Tk ■ If we now return to the quasiparticle density of states p* , we find it has 
the value 



at the Fermi energy so the mass enhancement of the heavy electrons is then 

m* _ Q qD 



1 + 



m 



(98) 



(99) 



3. The charge of the f-electron. 



How does the f-electron acquire its charge? We have emphasized from the beginning that the 
charge degrees of freedom of the original f-electrons are irrelevant, indeed, absent from the physics 
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of the Kondo lattice. So how are charged f-electrons constructed out of the states of the Kondo 
lattice, and how do they end up coupling to the electromagnetic field? 

The large theory provides an intriguing answer. The passage from the original Hamiltonian 
73\\ to the mean-field Hamiltonian ([83j) is equivalent to the substitution 



^S^pU)c]pC,o. Vflcja + Vc]j,^. (100) 

In other words, the composite combination of spin and conduction electron are contracted into a 
single fermi field 

n 

J J 

The amplitude \' ~ — /j^t'-j, = — ~ (t-'^jj /jj) involves electron states that extend over decades of 
energy out to the band edges. In this way the f-electron emerges as a composite bound-state of a 
spin and an electron. More precisely, in the long-time correlation functions. 

Such "clustering " of composite operators i nto a s i ngle e ntity is well known statistical mechanics 



as part of the operator product expansion (ICardvl . 



19961 ). In many body physics, we are used to 



the clu sterin g of fermions pairs into a composite boson, as in the BCS model of superconductiv- 
ity, —gy'^{x)i!>i{x') ^ A{x - x') . The unfamiliar aspect of the Kondo effect, is the appearance of a 
composite fermion. 

The formation of these composite objects profoundly modifies the conductivity and plasma 
oscillations of the electron fluid. The Read-Newns path integral has two U{1) gauge invariances 
- an external electromagnetic gauge invariance associated with the conservation of charge and 
an internal gauge invariance associated with the local constraints. The f-electron couples to the 
internal gauge fields rather than the external electromagnetic fields, so why is it charged? 

The answer lies in the broken symmetry associated with the development of the amplitude V. 
The phase of V transforms under both internal and external gauge groups. When V develops 
an amplitude, its phase does not actually order, but it does develop a stiffness which is sufficient 
to lock the internal and external gauge fields together so that at low frequencies, they become 
synonymous. Written in a schematic long-wavelength form, the gauge-sensitive part of the Kondo 
lattice Lagrangian is 

£ = Y, fd^'x [cUx)i-idt + e<^{x) + e^_^^)c„{x) + fl{x){-idt + A(x))^(x) 
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+ (V{x)cUx)Ux)+R.c. 



(103) 



where p = — iV. Suppose V{x) = a{x)e^'^^^\ There are two independent gauge transformations 
that that increment the phase (j) of the hybridization. In the external, electromagnetic gauge 
transformation, the change in phase is absorbed onto the conduction electron and electromagnetic 
field, so if y ^ Fe^", 

(j)^4> + (^, c(x) ^ c(x)e~*"(^\ e$(x) ^ e$(x) + q(x), eA ^ eA-Va{x). 

where ($, A) denotes the electromagnetic scalar and vector potential at site j and a = dta = —idrOi 
denotes the derivative with respect to real time t. By contrast, in the internal gau ge transformation 



the ph ase change of V is absorbed onto the f-fermion and the internal gauge field ([Read and Newns 



1983a|), so if y ^ Ve'l^, 

+ /(a;) ^ /(x)e^^(-), A(x) ^ A(x) - (104) 

If we expand the mean-field Free energy to quadratic order in small, slowly varying changes in 
A(x), then the change in the action is given by 



SS = -^ j d^xdT6X{x) 



where XQ = —6'^F/5X'^ is the f-electron susceptibility evaluated in the mean-field theory. However, 
6X{x) is not gauge invariant, so there must be additional terms. To guarantee gauge invariance 
under both the internal and external transformation, we must replace SX by the covariant combi- 
nation 5X + (j) — e^. The first two terms are required for invariance under the internal gauge group, 
while the last two terms are required for gauge invariance under the external gauge group. The 
expansion of the action to quadratic order in the gauge fields must therefore have the form 



S ~ J dT ^(0 + 5X{x) - e$(x))^ 



((jj -r u/\(X) — a^yx)^"^ 

j 

SO the phase (j) acquires a rigidity in time that generates a "mass" or energy cost associated with 
difference of the external and internal potentials. The minimum energy static configuration is 
when 

a(x)+0(x) = e«>(x), 

so when the external potential changes slowly, the internal potential will track it. It is this effect 
that keeps the Kondo resonance pinned at the Fermi surface. We can always choose the gauge 

47 




A(x) 




FIG. 19 (a) Spin liquid, or local moment phase, internal field A decoupled from electromagnetic field (b) 
Heavy electron phase, internal gauge field "locked" together with electromagnetic field. Heavy electrons are 
now charged and difference field [eyl(a:;) — A{x)] is excluded from the material. 



where the phase velocit y cj) is absorbed into t he local gauge field A. Recent work by Coleman, 



Marston and Schofield (Coleman et al 



2005al ) has extended this kind of reasoning to the case 
where RKKY couplings generate a dispersion jp-_A_ for the spinous, where A is an internal vector 
potential which suppresses currents of the gauge charge Uf. In this case, the long- wavelength action 
has the form 



S = ^ I cfxdT 



Ps[eA + V(P-AY - xqie^ - 4> - ^>^f 



In this general form, heavy electron physics can be seen to involve a kind of "Meissner effect" that 
excludes the difference field eA — A from within the metal, locking the internal field to the external 
electromagnetic field, so that the f-electrons which couple to it, now become charged (Fig. I19|) . 



4- Optical Conductivity of the heavy electron fluid. 



One of the interesting interesting consequences of the heavy electr on charge, is a complete 



renormalization of the electronic plasma frequency ((Millis et al 



1987bl ). The electronic plasma 



frequency is related via a f-sum rule, to the integrated optical conductivity 







— o-(w) = /l = - 
vr z 



Tire 



m 



(105) 
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where rig is the density of electrons. ^ In the absence of local moments, this is the total spectral 
weight inside the Drude peak of the optical conductivity. 

When the the heavy electron fluid forms, we need to consider the plasma oscillations of the 
enlarged Fermi surface. If the original conduction sea was less than half filled, then the renormalized 
heavy electron band is more than half filled, forming a partially filled hole band. The density of 
electrons in a filled band is N/ , so the effective density of hole carriers is then 

nHF = {N-Q- K)/a^ = (N - Q)/a^ - n,. 

The mass of the excitations is also renormalized, m — > m* . The two effects produce a low frequency 
'quasiparticle" Drude peak in the conductivity, with a small total weight 

'^^ ^ r \ f TTUHFe^ „ rn ( nHF\ „ ..„„x 

dua{u) = f2 = /i X — « /i (106 

2 m* m* \ Tic J 







Optical conductivity probes the plasma excitations of the electron fluid at low momenta. The direct 
gap between the upper and lower bands of the Kondo lattice are separated by a direct hybridization 
gap of order 2V ~ DTk ■ This scale is substantially larger than the Kondo temperature, and it 
defines the separation between the thin Drude peak of the heavy electrons and the high-frequency 
contribution from the conduction sea. 

In other words, the total spectral weight is divided up into a small "heavy fermion" Drude peak, 
of total weight /2, where 

(^i^) = — —-np^ — - (10^) 

m* (r* j ^ — luj 

separated off by an energy of order V ~ \/TkD from a n "inter-band" component associa t ed wit h 



excitations between the lower and upper Kondo bands (jDegiorgi et al 



2001 



Millis et al 



1987al ). 



This second term carries the bulk ~ /i of the spectral weight (Fig. [501 " 



Si mple calculations, b ased on the Kubo formula confirm this basic expectation, (iDegiorgi et al. 



2001 



Millis et al 



1987al ) showing that the relationship between the original relaxation rate of the 



The f-sum rule is a statement about the instantaneous, or short-time diamagnetic response of the metaL At short 
times dj/dt — {rice^ /■m)E, so the high frequency hmit of the conductivity is <j{lu) — ^^jr^- But using the 
Kramer's Kronig relation 

dx f^ix) 



a(cj) 

endequation at large frequencies, 



in X — uj — iS 



— iLO J n 

so that the short-time diamagnetic response implies the f-sum rule. 
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O((0) 




FIG. 20 Separation of the optical sum rule in a heavy fermion system into a high energy "inter-band" 
component of weight /2 ~ ne^ /m and a low energy Drude peak of weight /i ~ ne^ /m* . 



conduction sea and the heavy electron relaxation rate r* is 



m 



(108) 



Notice that this means that the residual resistivity 



m 



m 



Po 



(109) 



is unaffected by the effects of mass renormalization. This can be understood by observing that 
the heavy electron Fermi velocity is also renormalized by the effective mass, v*p = so that the 
mean- free path of the heavy electron quasiparticles is unaffected by the Kondo effect. 



I* = V*fT* 



VpT. 



(110) 



The formation of a narrow Drude peak, and the presence of a direct hybri d ization gap, have been 



seen in optical measure i nents on heavy electron systems (iBeyerman et al 



2001 



Schlessinger et al. 



19881 : 



Dordevic et al. 



19971 ). One of the interesting features about the hybridization gap of size 



2V, is that the mean-field theory predicts that the ratio of the direct, to the indirect hybridization 
gap is given by ~ "T^Y ~ ^Z^' ^° effective mass of the heavy electrons should scale 

as square of the ratio between the hybridization gap and the characteristic scale T* of the heavy 
Fermi liquid 



m 



■nif, 



oc 



2y Y 



50 



1,000 



100 



10 - 



0.1 




®CeAI, 

Frequency (cm ) ' 



0.1 



1 



10 

m /m. 



100 



1.000 



FIG. 21 Scaling of the effe ctive mass of heavy electrons with the square of the optical hybridization gap 



after (Dordevic et 



200 ih 



In practical experiments, Tk is replaced by the "coherence temperature" T* where the resistivity 
reaches a maxirnum. This scaling law is broadly followed (see Fig. [2T]) in measured optical data 



( Dordevic et al. 



2001 



), and provides provides further confirmation of the correctness of the Kondo 



lattice picture. 



D. Dynamical Mean Field Theory. 

The fermionic large N approach to the Kondo lattice provides a invaluable description of heavy 
fermion physics, one that can be improved upon beyond the mean-field level. For example, the 
fluctuations around the mean-field theory can be used to cq i npute the interactions, the dynamical 



correlation functions and the optical conductivity (jColemanl . 



1987b 



Millis et al. 



1987al ). However, 



the method does face a number of serious outstanding drawbacks: 

• False phase transition. In the large limit, the cross-over between the heavy Fermi liquid 
and the local moment physics sharpens into a phase transition where the 1/A^ expansion 
becomes singular. There is no known way of eliminating this feature in the 1/A^ expansion. 

• Absence of magnetism and superconductivity. The large approach, based on the SU (N) 
group, can not form a two-particle singlet for N > 2. The SU{N) group is fine for particle 
physics, where baryons are bound-states of N quarks, but for condensed matter physics, we 
sacrifice the possibility of forming two-particle or two-spin singlets, such as Cooper pairs and 
spin-singlets. Antiferromagnetism and superconductivity are consequently absent from the 
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mean-field theory. 



Amongst the various alternative approaches currently under consideration, one of particular 
note is the dynamical mean- field theory (DMFT). The idea of dynamical mean-field theory (DMFT) 



is to reduce the lattice problem to the p hysics of a sing ^ 



consistently determined effective medium (jGeorges et al. 



e magnetic ion embedded within a self- 



1996 



Kotliar et al 



20061). The effective 



medium is determined self-consistently from the self-energies of the electrons that scatter off the 
single impurity. In its more advanced form, the single impurity is replaced by a cluster of magnetic 
ions. 

Early versions of the dy namical mean-fiel d theory were cons idered by Kuromoto 



(jKuromoto and Watanabd . 



1971 



), Grewe and Cox (jCox and Grewd . Il988l ) and others, who used 



diagrammatic means to extract the physics of a single im purity. The modern conceptu al frame- 



work for DM FT was developed by Metzn er and Vollhardt (jMetzner and Vollhardt 



and Georges ([Georges and Kotliai 



Luttinger and Ward (IKotliar et al. 



19891 ') ■ Kothar 



19921) The basic idea behind DM FT is linked to early work of 



2006 



Luttinger and Ward . 



1960), who found a way of writing 



the Free energy as a variational functional of the full electronic Green's function 



(111) 



Luttinger and Ward showed that the Free energy is a variational functional of F[Q] from which 
Dyson's equation relating the Q to the bare Green's function Qq 

The quantity Ti[Q] is a functional, a machine which takes the full propagator of the electron and 
outputs the self-energy of the electron. Formally, this functional is the sum of the one-particle 
irreducible Feynman diagrams for the self-energy: while its output depends on the input Greens 
function, the actual the machinery of the functional is determined solely by the interactions. The 
only problem is, that we don't know how to calculate it. 

Dynamical mean-field theory solves this problem by approximating this functional by that of a 
single impurity or a cluster of magnetic impurities ( Fig. [22]). This is an ideal approximation for 



a local Fermi liquid, where the physics is highly retarded in time, but loca' 



proxi mation is also asymptotically exact in the limit of infinite dimensions ([Metzner and Vollhardt! . 



in space. The local ap' 



19891 ). If one approximates the input Green function to S by its on-site component Qij ~ Q^ij, 



then the functional becomes the local self-energy functional of a single magnetic impurity, 



(112) 
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impurity 

FIG. 22 In the Dynamical Mean Field Theory, the many body physics of the lattice is approximated by a 
single impurity in a self-consistently determined environment. Each time the electron makes a sortie from 
the impurity, its propagation through the environment is described by a self-consistently determined local 
propagator Giuj), represented by the thick red line. 



Dynamical mean field theory extracts the local self-energy by solving an Anderson impurity 
model embedded in an arbitrary electronic environment. The physics of such a model is described 
by a path integral with the action 

S = - dTdr'fUr)goHr-r')UT') + U f dTn^{T)n^{T) 
Jo Jo 

where Gq{t) describes the bare Green's function of the f-electron, hybridized with its dynamic 

environment. This quantity is self-consistently updated by the DMFT. There are by now, a large 

number of superb numerical methods to solve an Anderson model for an arbitrary environment, 

including the use of exac t diago nalization, diagrammatic techniques and the use of Wilson's renor- 



malization group (jBullal . 



20061 ) . Each of these methods is able to take an input "environment" 
Green's function providing as output, the impurity self-energy T,[Qq] = 

Briefly, here's how the DMFT computational cycle works. One starts with an estimate for the 
environment Green's function Qq and uses this as input to the "impurity solver" to compute the 
first estimate S(ia;„) of the local self-energy. The interaction strength is set within the impurity 
solver. This local self-energy is used to compute the Green's functions of the electrons in the 
environment. In an Anderson lattice, the Green's function becomes 



uj-E 



f 



w - ek 



-1 -1 



(113) 
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FIG. 23 Resist ivity for the Anderson lattice, calculated using the DMFT, computing the self-energy to 



order [/^, after ([Schweitzer and Czvchol] 



19911 ). (1), (2), (3) and (4) correspond to a sequence of decreasing 



electron density corresponding to utot = (0.8, 0.6, 0.4, 0.2) respectively. 

where V is the hybridization and ek the dispersion of the conduction electrons. It is through 
this relationship that the physics of the lattice is fed into the problem. From G(k, w) the local 
propagator is computed 



H -1 



io-E 



f 



w - ek 

Finally, the new estimate for the bare environment Green's function Qq, is then obtained by in- 
verting the equation = Gq^ — ^, so that 



(114) 



This quantity is then re-used as the input to an "impurity solver" to compute the next estimate 
of ^,{0;). The whole procedure is t hen r e-iterated to self-consistency. For the Anderson lattice, 



Cyzcholl (jSchweitzer and Czvcholll . 



1991 



y) has shown that remarkably good results are obtained 
using a perturbative expansion for S to order (Fig. [23l) . Although this approach is not 
sufficient to capture the limiting Kondo behavior, much the qualitative physics of the Kondo 
lattice, including the development of coherence at low temperatures is already captured by this 
approach. However, to go to the strongly correlated regime where the ratio of the interaction to the 
impurity hybridization width U /(ttA) is much larger than unity, one requires a more sophisticated 
solver. 

There are many ongoing developments under way using this powerful new computational tool, 
including the incorporation of realistic descriptions of complex atoms, and the extension to "cluster 
DMFT" involving clusters of magnetic moments embedded in a self-consistent environment. Let 
me end this brief summary with a list of a few unsolved issues with this technique 
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FIG. 24 Temperature dependent susceptibility in FeSi after paccarino et a/.l . 119671 ). fitted to the activated 
Curie form x(r) = {C/T)e-'^/'^''^^\ with C = igfiB)^jij + 1), and g ^ 3.92, A = 750K. The Curie tail 
has been subtracted. 



• There is at present, no way to relate the thermodynamics of the bulk to the impurity ther- 
modynamics. 

• At present, there is no natural extension of these methods to the infinite U Anderson or 
Kondo models that incorporates the Green's functions of the localized f-electron degrees of 
freedom as an integral part of the dynamical mean-field theory. 

• The method is largely a numerical black box, making it difficult to compute microscopic 
quantities beyond the electron-spectral functions. At the human level it is difficult for 
students and researchers to separate themselves from the ardors of coding the impurity 
solvers, and make time to develop new conceptual and qualitative understanding of the 
physics. 

III. KONDO INSULATORS 
A. Renormalized Silicon 



The ability of a dense lattice of local moments to transform a metal into an insulator, a 
"Kondo insulator" is one of the remarkable and striking consequences of the dense Kondo effect 
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FIG. 25 Schematic band pi cture of Kondo insulato r, illustrating how a magnetic field drives a metal insulator 
transition. Modified from (jAeppli and Fiskl . Il992[ l. 



(jAeppli and Fisk . 



1992 



Riseborough 



200C 



Tsunetsugu et al. 



19971 ). Kondo insulators are heavy 



electron systems in which the the liberation of mobile charge through the Kondo effect gives rise to 
a filled heavy electron band in which the chemical potential lies in the middle of the hybridization 
gap. From a quasiparticle perspective, Kondo insulators are highly renormalized "band insulators" 
(Fig. [25]) . The d-electron Kondo insulator FeSi has been referred to as "renormalized silicon". 
However, like Mott-Hubbard insulators, the gap in their spectrum is driven by interaction effects, 
and they display optical and magnetic properties that can not be understood with band theory. 
There are about a dozen known Ko n do iii sulators, inclu ding the rare earth s ystems SuiBq 



Menth et al 



( Hundlev et al. 



19691'! ■ YBi2 (llga et al 



Takabatake et al 



19901) . CeNiSn (jlzawa et al 



1988) CeFe^PT) (IMeisner a/, 



199S 



Takabatake et al 



1992 



1985|), CesBuPts 



1990) and CeRhSb 



1994 ). and the d-electron Kondo insulator FeSi. At high temperatures, Kondo 



insulators are local moment metals, with classic Curie susceptibilities, but at low temperatures, 
as the Kondo effect develops coherence, the conductivity and the magnetic susceptibility drop to- 
wards zero. Perfect insulating behavior is however, rarely observed due to difficulty of eliminating 
impurity band formation in ultra-narrow gap systems. One of the cleanest e xamples of Kondo insu - 



lating behavior occurs in the d-electron system FeSi (jPiTusa et al 



19971 : 



Jaccarino et al 



19671) • 



This "fly-weight" heavy electron system provides a rather clean realizati on of the phenomena see n 



in other Kondo insulators, with a spin and charge gap of about 750K (ISchlessinger et al 



19971 ). 



Unlike its rare-earth counterparts, the small spin-orbit coupling in this materials eliminates the 



Van Vleck contribution to the susceptibl 



completely vanishes at low temperatures (jJaccarino et al 



ity at T = 0, g i ving 



rise a susceptibility which almost 



19671 ) (Fig. (IMl). 



Kondo insulators can be understood as "half-filled" Kondo lattices in which each quenched 
moment liberates a negatively charged heavy electron, endowing each magnetic ion an extra unit 
of positive charge. There are three good pieces of support for this theoretical picture: 
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Each Kondo insulator has its fully itinerant semiconducting analog. For example, CeNiSn is 
iso-structural and iso-electronic with the semiconductor TiNiSi containing Ti^'^ ions, even 
though the former contains Ce'^"^ ions with localized / moments. Similarly, CesBi^Pt^, with 
a tiny gap of order lOmeV is isolectronic with semiconducting Th^Sh/i^Ni^, with a 70meV 
gap, in which the 5/ -electrons of the Th'^^ ion are entirely delocalized. 

Replacing the magnetic site with iso-electronic non-magnetic ions is equivalent to doping, 
thus in Cei-xLaxBi^Pts, each La^^ ion, behaves as an electron donor in a lattice of effective 
Ce^+ ions. Ces-xLaxPt^Bi^ is in fact, very similar to CePd^ which co ntains a pseudo-gap 



in its optical conductivity, with a tiny Drude peak (jBucher et al. 



1993 ) 



• The magneto-resistance of Kondo insulators is large and negative, and the "insulating gap" 
can be closed by the action of physically accessible fields. Thus in CesBi^Pts a 30T field is 
sufficient to close the indirect hybridization gap. 

These equivalences support the picture of the Kondo effect liberating a composite fermion. 

Fig. (|26l (a) ) shows the sharp rise in the resistivity of Ce^iBi^Pt^, as the Kondo insulating 
gap forms. In Kondo insulators, the complete elimination of carriers at low temperatures is also 
manifested in the optical conductivity. Fig. ([261 (b)) shows the temperature dependence of the 
optical conductivity in CesBi^Pt^, showing the emergence of a gap in the low frequency optical 
response and the loss of carriers at low energies. 

The optical conductivity of the Kondo insulators is of particular interest. Like the heavy 
electron metals, the development of coherence is marked by the formation of a direct hybridization 
gap in the optical conductivity. As this forms, a pseudo-gap develops in the optical conductivity. 
In a non-interacting band gap, the lost f-sum weight inside the pseudo-gap would be deposited 
above the gap. In heavy fermion metals, a small fraction of this weight is deposited in the Drude 
peak - however, most of the weight is sent off to energies comparable with the band-width of the 
conduction band. This is one of the most direct pieces of evidence that the formation of Kondo 
singlets involves electron energies that spread out to the bandwidth. Another fascinating feature of 
the heavy electron "pseudo-gap" , is that it forms at a temperature that is significantly lower than 
the pseudogap. This is because the pseudogap has a larger width given by the geometric mean of 
the coherence temperature and the band-width. 2V ~ \/TkD. The extreme upward transfer of 
spectral weight in Kondo insulators has not yet been duplicated in detailed theoretical models. For 
example, while calculations of the optical conductivity within the dynamical mean field theory, do 
show spectral weight transfer, it is not yet possible to reduce the indirect band-gap to the point 
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FIG. 26 (a) Tempe rature dependent res istivity of Ce-iPtiBi^ showing the sharp rise in resistivity at low 
temperatures after (jHundlev et all 11990 ) . Replacement of local moments with spinless La i ons acts like a 



dopa nt, (b) Real part of optical rconductivity <ti{lj) for Kondo insulator Ce^PtiBi^aliev (jBucher et al 



19941 ). The formation of the pseudogap associated with the direct hybridization gap, leads to the transfer of 
f-sum spectral weight to high energies of order the band-width. The pseudogap forms at temperatures that 
are much smaller than its width (see text). Insert shows (7i{uj) in the optical range. 



where it is radically smaller than the interaction scale U (iRozenberg et a/.l . Il996l ). It may be that 
these discrepancies will disappear in future calculations based on the more extreme physics of the 
Kondo model, but these calculations have yet to be carried out. 

There are however, a number of aspects of Kondo insulators that are poorly understood from 
the the simple hybridization picture, in particular 

• The apparent disappearance of RKKY magnetic interactions at low temperatures. 



The nodal character of the hybridization gap that develops in the narrowest gap Kondo 
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insulators, CeNiSn and CeRhSb. 



The nature of the metal-insulator transition induced by doping. 



B. Vanishing of RKKY interactions 



There are a number of experimental indications that the low-energy magnetic interactions vanish 
at low frequencies in a K ondo lattice. The low temperature product of the susceptibility and 



temperature xT reported ( Aeppli and Fis^ 



3, 



I992I ) to scale with the inverse Hall constant 1/Rh, 



representing the exponentially suppressed density of carriers, so that 



1 



-A/T 



X 



RhT T 

The important point here, is that the activated part of the susceptibility has a vanishing Curie- Weiss 
temperature. A similar conclusion is reached from inelastic neutron scattering measurements of the 
magnetic susceptibility x' {q,uj) ~ in CeNiSn and FeSi, which appears to lose all of its momentum 
dependence at low temperatures and frequencies. There is to date, no theory that can account for 
these vanishing interactions. 



C. Nodal Kondo Insulators 

The narrowest gap Kondo insulators, CeNiSn and CeRhSb are effectively semi-metals, for 
although they do display tiny pseudogaps in their spin and c harge spectra, the purest samples of 



these materials develop metallic behavior (|lzawa et al. 



19991 ). What is particularly peculiar (Fig. 
27|l about these two materials, is that the NMR relaxation rate l/(Ti)shows a temperature 



dependence from about lOK to IK, followed by a linear Korringa behavior at lower temperatures. 
The usual rule of thumb, is that the NMR relaxation rate is proportional to a product of the 
temperature and the thermal average of the electronic density of states A^*(a;) 

1 



Ti 



TN{loY ~ T[N{uj ~ T)Y 



(115) 



^I^^ ^ N{ujY is the thermally smeared average of the squared density of 



where N{ujf = f de y — g;^ 
states. This suggests that the electronic density of states in these material s has a "V shaped form, 
with a finite value at uj 



0. Ikeda and Miyake (jlkeda and Mivakd . 



19961 ) have proposed that the 



Kondo insulating state in these materials develops in a crystal field state with an axially symmetric 
hybridization vanishing along a single crystal axis. In such a picture, the finite density of states 
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does not derive from a Fermi surface, but from the angular average of the coherence peaks in the 
density of states. The odd thing about this proposal, is that CeNiSn and CeRhSb are monoclinic 
structures, and the low- lying Kramers doublet of the f-state can be any combination of the | ± 
I ± |) or I ± |) states: 

|± = 6i| ± 1/2) + 62I ± 5/2) + 63I ^ 3/2) 

where where b = (61, ^2 1^3) could in principle, point anywhere on the unit sphere, depending on 
details of the monoclinic crystal field. The Ikeda Miyake model corresponds to three symmetry- 
related points in the space of crystal field ground states. 



V2 _V5 3> 



4 ' 4 ' 4^ 

(0,0,1) 



(116) 



where a node develops along the x, y or z axis respectively. But the nodal crystal field states 
are isolated "points" amidst a continuum of fully gapped crystal field states. Equally strangely, 
neutron scattering results show no crystal field satellites in the dynamical spin susce ptibility of 



CeNiSn, suggesting that the crystal electr ic fields are quenched (lAlekseev et al 



the nodal physics is a many body effect (jKagan et al. 



19931 : 



Moreno and Coleman 



1994') ■ an d that 



2000). One 



idea, is that H und's interactions prov ide the driving force for this selection mechanism. Fulde 



and Zwicknagl (jZwicknagl et al 



2002l )have suggested that Hund's coupling 's select the orbitals in 



multi f-electron heavy electron metals such as UPt^. Moreno and Coleman (jMoreno and Colemanl . 



2OOOI ) propose a similar idea, in which virtual valence fluctuations into the state generate a 
many-body, or Weiss effective field that couples to the orbital degrees of freedom, producing an 
effective crystal field which adjusts itself in order to minimize the kinetic energy of the f-electrons. 
This hypothesis is consistent with the observation that the Ikeda Miyake state corresponds to the 
Kondo insulating state with the lowest kinetic energy, providing a rational for the selection of the 
nodal configurations. Moreno and Coleman also found another nodal state with a more marked 
^-shaped density of states that may fit the observed properties more precisely. This state is also a 
local minimum of the electron Kinetic energy. These ideas are still in their infancy, and more work 
needs to be done to examine the controversial idea of a Weiss crystal field, both in the insulators 
and in the metals. 
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FIG. 27 (a) NMR relaxation rate 1/Ti in CeRhSb and CeNiSn, showing a relaxation rate sandwiched 
between a l ow, and a high tempera ture T-linear Korringa relaxation rate, suggesting a V- shaped density of 



states after ( Nakamura et al. 



1994 ) (b) Contour plot of th e grou nd-state energy in mean-field theory for the 



narrow gap Kondo insulators after (jMoreno and Colemanl . uOOOT) , as a function of the two first components 
of the unit vector b (the third one is taken as positive). The darkest regions correspond to lowest values of 
the free-energy. Arrows point to the three global and three local rninima that correspond to nodal Kondo 
insulators. (c) Density of states of Ikeda Miyake (jlkeda and Mivakel Il99(: ) state that appears as the g lobal 
minimum of the Kinetic energy (d) Density of states of the MC state (jMoreno and Colemanl . |2000() that 
appears as a local minimum of the Kinetic energy, with more pronounced "V" shaped density of states. 

IV. HEAVY FERMION SUPERCONDUCTIVITY 



A. A quick tour 



Since the discovery (jSteglich et al. 



1970 ) of superconductivity in CeCu2Si2, the list of known 



heavy fermion superconductors has grow n to include more a doz en 



1995 



Sigrist and Ueda 



Sigrist and Ueda 



1991b|) mate- 



1991al ). In each 



rials with a great diversity of properties (jCox and Maple 
of these materials, the jump in the specific heat capacity at the superconducting phase transition 
is comparable with the normal state specific heat 



Cv 



1 



(7) 



and the integrated entropy beneath the C^/T curve of the superconductor matches well with the 
corresponding area for the normal phase obtained when superconductivity is suppressed by disorder 
or fields 



I 

Jo 



dT- 



0. 



(8) 
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Since the normal state entropy is derived from the f-moments, it foHows that these same degrees 
of freedom are involved in the development of the superconducting state. With the exception of a 
few anomalous cases, {UBei^, PuCoGa^ and CeCoIn^) heavy fermion superconductivity develops 
out of the coherent, paramagnetic heavy Fermi liquid, so heavy fermion superconductivity can be 
said to involve the pairing of heavy f-electrons. 

Independent confirmation of the "heavy" nature of the pairing electrons comes from observed 
size of the London penetration depth and superconducting coherence length ^ in these systems, 
both of which reflect the enhanced effective mass. The large mass renormalization enhances the 
penetration depth, whilst severely contracting the coherence length, making these extreme type-II 
superconductors. The London penetration depth of heavy fermion superconductors agree well with 
the value expected on the assumption that only spectral weight in the quasiparticle Drude peak 
condenses to form a superconductor by 

1 ne^ f doj , . ne^ 

-a{u)) « (9) 



London penetration depths in these compounds are a factor of 20 — 30 times longer (iBroholm et al. 



19901 ) than in superconductors, reflecting the large enhancement in effective mass. By contrast, 
the coherence lengths vp/ A ^ hkp / {m* A) are severely contracted in a heavy fermion super- 
conductor. The orbitally limited upper critical field is determined by the condition that an area 

contains a flux quantum (,^Bc ~ In UBei^, a superconductor with 0.9K transition temper- 
ature, the upper critical field is about 11 Tesla, a value about 20 times larger than a conventional 
superconductor of the same transition temperature. 

Table B. shows a selected list of heavy fermion superconductors. "Canonical" heavy fermion 
superconductors, such as CeCu2Si2 and UPt^, develop superconductivity out of a paramag 



netic Landau 



Geibel. et al 



erm i liquid. "Pre-ordered" superconductors, such as U Pt2Al^ (jGeibel et al. 



1991 



1991 



, CePt-^Si and URu2Si2, develop another kind of order before going supercon- 
ducting at a lower temperature. In the case of URu2Si2, the nature of the upper ordering transition 
is still unidentified, but in the other examples, the upper transition involves the development of an- 
tiferromagnetism. "Quantum critical" superconductors, including Celn^ and CeCu2{Sii-xGex)2 
develop superconductivity when pressure tuned close to a quantum critical point. Celn^ devel- 
ops superconductivity at the pressure-tuned antiferromagnetic quantum critical point at 2.5GPa 
{25kbar). GeCu2{Si, Ge)2 has two islands, one associated with antiferromagnetism at low pressure, 
and a secon d at s till higher pressure, thought to be associated with critical valence fluctuations 



(jYuan et al. 



20061). 
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Table. B. Selected Heavy Fermion Superconductors. 



Type 




T 

J- c 


Knio-ht mi iff 
(Singlet) 




Can 
Symmetry 




"Canonical" 


CeCu2Si2 


0.7 


Singlet 


First HFSC 


Line nodes 


[1] 




0.48K 


? 


Double Transition to 
T-violating state 


Line and 
point nodes 


[2] 


Pre-ordered 


UPd2Ah 


2K 


Singlet 


IN eel Ai^iVl 
Tn = UK 


Line nodes 
A ~ cos 2x 


[3] 


U Ru2Sl2 


1.3K 


Singlet 


Hidden order at 
To = 17.5K 


Line nodes 


[4] 


CePtsSi 


0.8K 


Triplet 


\-^^^T'^i'\^—'\^^^\\ si I'm c "vl" si 1 
± diiLV viUiajLiiig Ji-Lcxi. 

Tn = 3.7 K 


Line nodes 


[5] 


Quantum 
Critical 


Celn^ 


(2.5GPa) 


Singlet 


±. 1± O t IJU.ClilLjU.ili l^ilLH^Cli 

HFSC 


Line nodes 


[6] 


CeCu2{Sii^xGex)2 


iK (P—0) 
0.95A'(5.4GPa) 


Singlet 


Two i^lanrl^ of SC^ ?m 
function of pressure 


Line nodes. 


[7] 


Quadrupolar 


PrOsiSbi2 


1.85K 


Singlet 


Quadrupolar 
fluctuations 


Point nodes 


[8] 


"Strange" 


CeCoIn^ 


2.3K 


Singlet 


(^uasi-zU 


Line nodes 


[9] 


UBei3 


OMK 


? 


Incoherent metal at Tc- 


Line nodes 


[10] 


PuCoGufj 


18.5K 


Singlet 


Direct transition Curie 
metal HFSC 


Line nodes 


[11] 




Steslid 


1 et al 


, W7( 




Stewart et al. 


. 1984R 


). f3] i Geibel et a 




1991: Sato et al.. 2001: 


Tou et al.. 




1991 


): [4] I 


Kim et al.. 


2003: 


Palstra e 


: al.. 


1985), [51 


Bauer et a 




2004), [6] 


iMathur et al.. 


199 




lYuan et al 




200 


d). 18] ilsawa et al.. 


200c 


). [9] ( 


Petrovic et al.. 


2001 


), [10] 1 


Andres et al.. 


1975 


, [11] 


i Sarrao et al. 


20 









"Strange" superconductors, which include UBeis, the 115 material CeCoIn^ and PuCoGa^, 
condense into the superconducting state out of an incoherent or strange metallic state. U Bei3 has a 
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resistance of order IbOfiflcm at its transition temperature. CeCoIn^ bears superficial resemblance 
to a high temperature superconductor, with a linear temperature resistance in its normal state, 
while its cousin, PuCoGa^ transitions directly from a Curie paramagnet of unquenched f-spins into 
an anisotropically paired, singlet superconductor. These particular materials severely challenge 
our theoretical understanding, for the heavy electron quasiparticles appear to form as part of the 
condensation process, and we are forced to address how the f-spin degrees of freedom incorporate 
into the superconducting parameter. 



B. Phenomenology 

The main body of theoretical work on heavy electron systems is driven by experiment, and 
focuses directly on the phenomenology of the superconducting state. For these purposes, it is 
generally sufficient to assume a Fermi liquid of pre- formed mobile heavy electrons, an electronic 
analog of superfluid Helium-3, in which the quasiparticles interact through a phenomenological BCS 
model. For most purposes. Landau Ginzburg theory is sufficient. I regret that in this short review, 
I do tno have time to properly represent and discuss the great wealth of fascinating phenomenology 
- the wealth of multiple phases, and the detailed models t hat have been developed to describe them. 
I refer the interested reader to reviews on this subject. (jSigrist and Uedal . Il991al ). 

On theoretical grounds, the strong Coulomb interactions of the f-electrons that lead to mo- 
ment formation in heavy fermion compounds are expected to heavily suppress the formation of 
conventional s-wave pairs in these systems. A large body of evidence favors the idea that the gap 
function and the anomalous Green function between paired heavy electrons F^pix) = (cq(x)c^(O)) 
is spatially anisotropic, forming either p-wave triplet, or d-wave singlet pairs. 

In BCS theory, the superconducting quasiparticle excitations are determined by a one-particle 
Hamiltonian of the form 



where 



Aa/3(k) = < 



A(k)(io-2)Q,/3 (singlet), 
d(k) • {ia2a)af3 (triplet). 



For singlet pairing, A(k) is an even parity function of k while for triplet pairing, (i(k) is a an 
odd-parity function of k with three components. 
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The excitation spectrum of an anisotropically paired singlet superconductor is given by 



This expression can also be used for a triplet superconductor that does not break time-reversal 
symmetry, by making the replacement |Akp = (ij^dk. 

Heavy electron superconductors are anisotropic superconductors, in which the gap function van- 
ishes at points, or more typically, along lines on the Fermi surface. Unlike s-wave superconductors, 
magnetic and non-magnetic impurities are equally effective at pair-breaking and suppressing Tc in 
these materials. A node in the gap is the result of sign changes in the underlying gap function. If 
the gap function vanishes along surfaces in momentum space, and the intersection of these surfaces 
with the Fermi surface produces "line nodes" of gapless qu asiparticle exc i tations. As an exam- 



ple, c onsider UPh, where, according to one set of models ([B 



1990 



Jovnt 



19881 : 



Machida and Ozakil . 



19891 : 



Puttika and Jovnt 



ount et al. 



199C 



Hess and Sauls 



19881 ). pairing involves a complex 



d-wave gap 



Here vanishes along the basal plane A;^ = 0, producing a line of nodes around the equator of 
the Fermi surface, and along the z-axis, producing a point node at the poles of the Fermi surface. 

One of the defining properties of line nodes on the Fermi surface is a quasiparticle density of 
states that vanishes linearly with energy 



N*{E) = 2^5{E - Ei,) ocE 



The quasiparticles surrounding the line node have a "relativistic" energy spectrum. In a plane 



perpendicular to the node, E]^ ~ {vpkiY + (0^2)^ where a = d/S./dk2 is the gradient of the gap 
function and ki and k2 the momentum measured in the plane perpendicular to the line node. For 
a two dimensional relativistic particle with dispersion E = ck, the density of states is given by 
N(E) = For the anisotropic case, we need to replace c by the geometric mean of vp and a, so 



vpa. This result must then be doubled to take account of the spin degeneracy and averaged 



over each line node: 



N{E) = 2 J] / 



dku 



\E\ 



27r ATTVpct 



1^1 X E 

nodes 



dki 



nodes ' 

In the presence of pair-breaking impurities, and in a vortex state the quasiparticle nodes are 
smeared, adding a small constant component to the density of states at low energies. 
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This linear density of states is manifested in a variety of power laws in the temperature depen- 
dence of experimental properties, most notably 

• Quadratic temperature dependence of specific heat Cy oc T^, since the specific heat coeffi- 
cient is proportional to the thermal average of the density of states 

ocT 

Cv 



^ oc N{E) ~r 



where N{E) denotes the thermal average of N{E). 

• A ubiquitous temperature dependence in the nuclear magnetic and quadrapole relaxation 
(NMR/NQR) relaxation rates 1/Ti. The nuclear relaxation rate is proportional to the 
thermal average of the squared density of states, so for a superconductor with line nodes, 

oc T N{EY ~ 

Fig. [28] shows the NMR relaxation rate of the Aluminum nucleus in UPd2Al^. 

Although power-laws can distinguish line and point nodes, they do not provide any detailed 
information about the triplet or singlet character of the order parameter, or the location of the 
nodes. The observation of upper critical fields that are "Pauli limited" (set by the spin coupling, 
rather than the diamagnetism) , and the observation of a Knight shift in most heavy fermion 
superconductors, indicates that they are anisotropically singlet paired. Three notable exceptions 
to this rule are UPt^, UBeis and UNi2Al^, which do not display either a Knight shift or a Pauli- 
limited upper critical field, are the best candidates for odd-parity triplet pairing. In the special 
case of CePt^Sn, the crystal structure lacks a center of symmetry and the resulting parity violation 
is must give a mixture of triplet and singlet pairs. 

Until comparatively recently, very little was known about the positions of the line nodes in heavy 
electron superconductors. In one exception, experiments carried out almost twenty years ago on 



UPt^ observed mar 



(jBishop et al 



i ed anisotropics i n the ultrasound attenuation length and the penetration depth. 



1984 



Broholm et al 



19901 ) that appear to support a line of nodes in the basal plane. 
The ultrasonic attenuation as{T)/an in single crystals of UPt^ has a T linear dependence when 
the polarization lies in the basal plane of the gap nodes but a dependence when the polarization 
is along the c-axis. 

An interesting advance in the experimental analysis of nodal gap structure has recently occurred, 
thanks to new insights into the behavior of the nodal excitation spectrum in the flux phase of heavy 
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FIG. 28 Temperature dependence of the ^'Al NQR relaxation rate l/Ti for UPd2Al^ after (jTou et a/. 



19951 ) showing T'^ dependence associated with hnes of nodes. I nset showing noda l structure A cx cos(2^) 



proposed from analysis of anisotropy of thermal conductivity in (|Won et al. 



2004 . 



fermion superconductors. In the nineties, Volovik (IVolovik 



19931 ) observed that the energy of heavy 



electron quasiparticles in a flux lattice is "Doppler shifted" by the superflow around the vortices, 
giving rise to a finite density of quasiparticle states around the gap nodes. The Doppler shift in 
the quasiparticle energy resulting from superflow is given by 



where Vg is the superfluid velocity and (p the superfluid phase. This has the effect of shifting 

'2R^ 



quasiparticle states by an energy of order AE' ~ where R is the average distance between 



^ is the flux 



vortices in the flux lattice. Writing t:HE? ~ <I>Oi and ttHc2S,'^ ~ where 
quantum, Hc2 is the upper critical field and ^ is the coherence length, it follows that ~ | 
Putting ^ ~ vp/ A, where A is the typical size of the gap, the typical shift in the energy of nodal 
quasiparticles is of order Eh ~ /S.^J'^^. Now since the density of states is of order N{E) = ■l^A^(O), 
where A^(0) is the density of states in the normal phase, it follows that the smearing of the nodal 
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quasiparticle energies will produce a density of states of order 



This effect, the "Volovik effect", produces a linear component to the specific heat Cy /T oc \J'^~- 
This enhancement of the density of states is largest when the group velocity Vp at the node is 
perpendicular to the applied field and when the field is parallel to iTf at a particular node, the 
node is unaffected by the vortex lattice (Fig. [29l) . This gives rise to an angular dependence in the 
specific heat coefficient and thermodynamics that can be used to measure the gap anisotropy. In 
practice, the situation is complicated at higher fields where the Andreev scattering of quasiparticles 
by vortices becomes important. The case of CeCoIn^ is of particular current interest. Analyses 
of the field-anisotropy of the thermal conductivity in this material was interpreted early on in 
terms of a gap structure with d^2_„2, while the anisotropy in the sp ecific heat appears to su g gest a 



20061) 



dxy symmetry. Recent theoretical work by Vorontsev and Vekhter (jVorontsov and Vekhteii . 
suggests that the discrepancy between the two interpretations can be resolved by taking into 
account the effects of the vortex quasiparticle scattering that were ignored in the specific heat 
interpretation. They predict that at lower fields, where vortex scattering effects are weaker, the 
sign of the anisotropic term in the specific heat will reverse, accounting for the discrepancy 

It is clear that, despite the teething problems in the interpretation of field- anisotropies in trans- 
port and thermodynamics, this is an important emerging tool for the analysis of gap anisotropy, 
and to date, it has been used to give tentative assignments to the gap anisotropy of UPd2Al^, 
CeCoIn^ and PrOsiShi2. 



C. Microscopic models 



1. Antiferromagnetic fluctuations as a pairing force. 



The classic theoretical models for heavy fermion superconductivity treat the heav y electron flu- 



ids as a Fermi liquid with antiferromagnetic interactions amon gst their quasiparticles (iMonod et al 



1986 



Monthoux and Lonzarichl . 



19991 : 



Scalapino et al 



19861 ). UPt^ provided the experimental in- 



spiration for early theories of heavy fermion superconductivity, for its superconducting state forms 
from within a well-developed Feri ni liquid. Neutron scattering on this material shows signs of an- 
tiferromagnetic spin fluctuations ((Aeppli et a/.l . 119871 ). making it natural to presuppose that these 
might be the driving force for heavy electron pairing. 
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v,*'- 'iJ "V^ 



CeCoIn g 

FIG. 29 Schematic sh owing how the nodal quasiparticle density of states depends on field orientation after 



(Vekhter et al 



19991 ). (a) Four nodes are activated when the field points towards an antinode, creating 



a maximum in density of states, (b) Two nodes activated when the field points towards a node, creat- 
i ng a minimum in th e density of states, (c) Theoretical dependence of density of s t ates o n angle, after 



( Vekhter et al 



19991) . (d) Measured angular dependence of C^/T after (lAoki et 



20041) is 45° out of 



phase with prediction. This discrepancy is believed to be due to vortex scattering, and is expected to vanish 
at lower fields. 

Since the early seventies, theoretical models had predicted that strong ferromagnetic spin fluc- 
tuations, often called "paramagnons" , could induce p-wave pairing, and it this mechanism was 
widely held to be the driving force for pairing in superfluid He — 3. An early proposal that anti- 
ferromagn etic interactio ns could provide the driving force for anisotropic singlet pairing was made 
19851). Short ly thereafter, three seminal papers. 



by Hirsch (iHirschl. 



and Emery (jMonod et al. 



jy Beal-Monod, Bour bonnais 



19861 ) (BEE), Sc alapino, Loh a n d Hir sch (|ScalaDino et al\ . Il986l ) (SLH) 



and by Miyake, Schmitt-Rink and Varma (jMivake et al 



1986|) (MSV), sohdified this idea with 



a concrete demonstration that antiferromagnetic interactions drive an attractive BCS interaction 
in the d-wave pairing channel. It is a fascinating thought that at the same time that this set of 
authors were forging the foundations of our current thoughts on the link between antiferromag- 
netism and d-wave superconductivity, Bednorz and Mueller were in the process of discovering high 
temperature superconductivity. 

The BBE and SLH papers develop a paramagnon theory for d-wave pairing in a Hubbard model 
with a contact interaction /, having in mind a system, which in the modern context, would be said 
to be close to an antiferromagnetic quantum critical point. The MSV paper starts with a model 
with a pre-existing antiferromagnetic interaction, which in the modern context would be associated 
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with the "t-J" model. It is this approach that I sketch here. The MSV model is written 

H = ^ ekOk^OkfT + Hint (117) 

where 

Hint = ^^^Ji^- kO-^a/? • ^-yS X (aL+q/2a"-k+q/27) ( «-k'+q/25ak'+q/2/3) (US) 
k,k' q 

describes the antiferromagnetic interactions. There are a number of interesting points to be made 
here: 

• The authors have in mind a strong coupled model, such as the Hubbard model at large U 
where the interaction can not be simply derived from paramagnon theory. In a weak-coupled 
Hubbard model a contact interaction I and bare suscept ibility xn(Q), the ind uced magnetic 



1989 1 as 



interaction can be calculated in an RPA approximation (|Mivake et al 

MSV make the point that the detailed mechanism that links the low-energy antiferromagnetic 
interactions to the microscopic interactions is poorly described by a weak-coupling theory, 
and quite likely to involve other processes, such as the RKKY interaction, and the Kondo 
effect that lie outside this treatment. 

Unlike phonons, magnetic interactions in heavy fermion systems can not generally be re- 
garded as retarded interactions, for they extend up to an energy scale loq that is comparable 
with the heavy electron band- width T* . In a classic BCS treatment, the electron energy 
are restricted to lie within a Debye energy of the Fermi energy. But here, ujq ~ T*, so all 
momenta are involved in magnetic interactions, and the interaction can transformed to real 
space as 

where J(R) = X]q e**''"^ J(q) is the Fourier transform of the interaction and ai = al^aa/BCiiff 
is the spin density at site i. Written in real space, the MSV model is seen to be an early 
predecessor of the t — J model used extensively in the context of high temperature super- 
conductivity. 
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To see that antiferromagnetic interactions favor d-wave pairing, one can use the, let us decouple 
the interaction in real space in terms of triplet and singlet pairs. Inserting the identity ^ 



3 1 



(120) 



into (jll9l ) gives 



H.. 



int 



(121) 



where 



(122) 



create singlet and triplet pairs with electrons located at sites i and j respectively. In real-space it 
is thus quite clear that an antiferromagnetic interaction Jij > induces attraction in the singlet 
channel, and repulsion in the triplet channel. Returning to momentum space, substitution of (jl20l 
) into (|118p gives 



H.nt = - 5^ J(k - k') qM/k', q - ^L, q • ^k', q 



k,k' 



(123) 



where ^-j^^ = i [a^^^/^ «(-^(^2)a7 «-k-q/2 y) ^nd = i „(-^<?^2) 
ate singlet and triplet pairs at momentum q respectively. Pair condensation is described by the 
zero momentum component of this interaction, which gives 



1 J 



1 J 



jay 0.-\i-c^/2 7 ' 



Hint = E [^kJ^Uk' + V^l,^i . ^k' 
k,k' 



(124) 



To prove this identity, first note that any two dimensional matrix, M can be expanded as M = moo'2 + rh ■ (J2(t, 
(6 = (1,3)) where mo = iTr[M(T2] and m = ^TT[Maa2], so that in index notation 

M^.f = iTr[Mcr2](cr2)„-, + iTr[A/<Tcr2] • (cr2a)c«7 
Now if we apply this relationship to the ay components of cfap ■ S-^s , we obtain 

b— 1,3 

If we now use the relation a'^a2 — ~a2ff, together with a • a = 3 and aabS = —at, we obtain 

- ^ 3 1 ^ 

Sap ■ S-yS ~ — -(0-2)0,7(0-2)^^ + -((TO-2)o7 • (o-2o)^^ 
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where = 5 (aL(~'^^2)a/3 ^-k/j) *k,q = 1 (aka(~'^'^'^2)a/3 O-k/?) create Cooper pairs and 



•^k.k' — 
yit) _ 



-3[J(k-k') + J(k + k')]/2 
[J(k-k')- J(k + k')]/2 



(125) 



are the BCS pairing potentials in the singlet, and triplet channel respectively. (Notice how the 
even/odd parity symmetry of the triplet pairs pulls out the corresponding symmetrization of J(k — 
k')) 

For a given choice of J(q), it now becomes possible to decouple the interaction in singlet and 
triplet channels. For example, on a cubic lattice of side length, if the magnetic interaction has the 
form 

J(q) = 2J{cos{qxa) + cos{qya) + cos(g^;a)) 

which generates soft antiferromagnetic fluctuations at the staggered Q vector Q = (7r/o,7r/a,7r/a), 
then the pairing interaction can be decoupled into singlet and triplet components, 

3J 



^k,k' = ^ [s(k)s(k') + 42_j,2(k)42_2^2(k') + d2^2_^2(k)d2z2_r2(k')] 

i=x,y,z 

where 

s(k) = y^(cos(A;x-a) + cos{kya) + cos(A;^a)) 
o?^2_^2(k) = (cos(A;^a) — cos{kya), 
(i222_^2(k) = ;^(cos(fe^a) + cos{kya) — 2cos{kza)) 

are the gap functions for singlet pairing and 



(126) 



(extended s-wave) 
(d-wave) 



(127) 



Pi{q) = V2sm{qia), {i = x, y, z), 



(p-wave) 



(128) 



describe three triplet gap functions. For J > 0, this particular BCS model will then give rise to 
extended s and d-wave superconductivity with approximately the same transition temperatures, 
given by the gap equation 

Vtanhfi^Ul , 
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2. Towards a unified theory of HFSC 



Although the spin fluctuation approach described provides a good starting point for the phe- 
nomenology of heavy fermion superconductivity, it leaves open a wide range of questions that 
suggest this problem is only partially solved: 

• How can we reconcile heavy fermion superconductivity with the local moment origins of the 
heavy electron quasiparticles? 

• How can the the incompressibility of the heavy electron fluid be incorporated into the theory? 
In particular, extended s-wave solutions are expected to produce a large singlet f-pairing 
amplitude, giving rise to a large Coulomb energy. Interactions are expected to significantly 
depress, if not totally eliminate such extended s-wave solutions. 

• Is there a controlled limit where a model of heavy electron superconductivity can be solved? 

• What about the "strange" heavy fermion superconductors UBei3, CeCoIn^ and PuCoGa^, 
where Tc is comparable with the Kondo temperature? In this case, the superconducting 
order parameter must involve the f-spin as a kind of "composite" order parameter. What is 
the nature of this order parameter, and what physics drives Tc so high that the Fermi liquid 
forms at much the same time as the superconductivity develops? 

One idea that may help to underst a,nd the heavy fe rmion pairing mechanism is Anderson's 



resonating valence bond (RVB) theorv (lAndersonl. 



Anderson proposed ([Anderson 



1987 



Baskaran et al 



19871) of high temper ature superconductivity. 



19871 : 



Kotliar 



19881 ) that the parent state of 



the high temperature superconductors is a two-dimensional spin liquid of resonating valence bonds 
between spins, w hich becomes superconductin g upon doping with holes. In the early nineties, Cole- 



man and Andrei (Coleman and Andrei 



19891 ) adapted this theory to a Kondo lattice. Although an 



RVB spin-liquid is unstable with respect to antiferromagnetic order in three dimensions, in situa- 
tions close to a magnetic instability, where the energy of the antiferromagnetic state is comparable 
with the Kondo temperature, Eafm ~ Tk, conduction electrons will partially spin-compensate 
the spin liquid, stabilizing it against magnetism (Fig. [30] (a)). In the Kondo-stabilized spin liquid, 
the Kondo effect induces some resonating valence bonds in the f-spin liquid to escape into the 
conduction fluid where they pair charged electrons to form a heavy electron superconductor. 
A key observation of RVB theory is that when charge fluctuations a re removed to form a spin 



fluid, there is no distinction between particle and hole ([Affleck et al. 



88l ). The mathematical 
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(a) 



Enirgy 



(c) 



F — T 

KoiHl4>-^t»l>lllzecl spin liquid .SY. K 



(b) I 





f— spins 



f— spi ns 



1.i> 




FIG. 30 Kondo stabilized spin liquid, diagram from (jColeman and Andreil . ll989l )(a)Spin liquid stabilized by 



Kondo effect, (b) Kondo effect causes singlet bonds to form between spin liquid and conduction sea. Escape 
of these bonds into the conduction sea induces superconductivity, (c) Phase diagram computed using SU{2) 
mean-field theory of Kondo Heisenberg model. 



consequence of this, is that the the spin-1/2 operator 



is not only invariant under a change of phase /o- — > e^'^fa, but it also possesses a continuous particle 
hole symmetry 



a 



fl cos Ofl + sgnasinOf^^. 



(129) 



These two symmetries combine to create a local SU (2) gauge symmetry. One of the implications 
is that the constraint nj = 1 associated with the spin operator, is actually part of a triplet of 
Gutzwiller constraints 



If we introduce the Nambu spinors 



f},h-fiifl=o, 44 = 0, MM = o. 



(130) 



fi 



fi ~ (4' fil)- 



(131) 



then this means that all three components of the "isospin" of the f-electrons vanish, 



1 

1 0, 



-i 

1 



1 

-1, 



(132) 
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where r is a triplet of Pauli spin operators that act on the f- Nambu spinors. In other words, in 
the incompressible f-fluid, there can be no s — wave singlet pairing. 

This symmetry is preserved in spin-1/2 Kondo models. When applied to the Heisenberg Kondo 
model 

ka («J') j 

where Si = fj^ ( S ) /j^ represents an f-spin at site i, it leads to an SU (2) gauge theory for the 
Kondo lattice with Hamiltonian 

1 

Jk 



where \j is the Lagrange multiplier that imposes the Gutzwiller constraint r = at each site, 



Cfc = ^J*"! J ^-nd Cj = (^ct^J ^'^^ Nambu conduction electron spinors in the momentum, and 
position basis, respectively, while 

are matrix order parameters associated with the Heisenberg and Kondo decoupling, respectively. 
This model has the local gauge invariance fj djfj, Vj djVj Uij — > giUijgj where gj is an SU{2) 
matrix. In this kind of model, one can "gauge fix" the model so that the Kondo effect occurs in 
the particle- hole channel (oj = 0). When one does so however, the spin-liquid instability takes 
place preferentially in an anisotropically paired Cooper channel. Moreover, the constraint on the 
f-electrons not only suppresses singlet s-wave pairing, it also suppresses extended s-wave pairing 
(Fig. [30]). 

One of the initial difficulties with both the RVB and the Kondo stabilized-spin liquid approaches 
is that in its original formulation, it could not be integrated into a large approach. Recent 
work indicates that both the fermionic RVB and the Kondo stabilized spin liquid picture can be 
formulat ed as a controlled SU(2) g auge theory by carrying out a large expansion using the group 







(V a^\ 






[ai -Vi) 



SP{N) (iRead and Sachdevl . Il99ll ). originally introduced by Read and Sachdev for problems in 
frustrated magnetism, in place of the group SU{N). The local particle-hole symmetry associated 
with the spin operators in SU{2) is intimately related to the symplectic property of Pauli spin 
operators 

(T2<?^o"2 = —a (135) 
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20021 ). The 



FIG. 31 Temperature dependence of the magnetic susccptibihty of PuCoGa^ after (jSarrao et al. 
susceptibiUty shows a direct transition from Curie Weiss paramagnet into heavy fermion superconductor, 
without any intermediate spin quenching. 

where is the transpose of the spin operator. This relation, which represents the sign reversal 
of spin operators under time-reversal, is only satisfied by a subset of the SU (N) spins for N > 2. 
This subset defines the generators of the symplectic subgroup of SU{N), called SP{N). 



Concluding this sectio n, I want to 



perconductor, PuCoGa^ (Curro et al 



jriefiy mention the chall enge posed by the highest Tc su- 



200 



1 



Sarrao et al 



2OO2I ). This material, discovered some 



four years ago at Los Alamos, undergoes a direct transition from a Curie paramagnet into a heavy 
electron superconductor at around Tc = 19K (Fig l3ip . The Curie pararnagneti sm is also seen 



in the Knight-shift, which scales with the bulk susceptibility (|Curro et al 



20051 ). The remark- 

able feature of this material, is that the specific heat anomaly has the large size [IWmJ/mol/ K'^ 



(Sarrao et al 



2002! )) characteristic of heavy fermion superconductivity, yet there there are no signs 



of saturation in the susceptibility as a precursor to superconductivity, suggesting that the heavy 
quasiparticles do not develop from local momen ts until the t ransit ion. This aspect of the physics 



can not be explained by spin-fluctuation theory llBang et al 



20041 ). and suggests that the Kondo 



effect takes place simultaneously with the pairing mechanism. One interesting possibility here, is 
that the development of coherence between the Kondo effect in two different channels created by 
the different symmetries of the valence fluctuations iri t o the and states rn ight be the driver 



of this intriguing new superconductor ([Coleman et al 



1999 



Jarrell et al 



19971) 
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V. QUANTUM CRITICALITY 



A. Singularity in the Phase diagram 

Many heavy electron systems can be tuned, with pressure, chemical doping or applied mag- 



Miranda and Dobrosavlievid . 



2005 



Stewart 



" ( 


Coleman et al. 




2001 




von Lohnevsen et al. 




2007; 


2001 




2006 




Varma et al 




2OO2I). The remarkable 



transformation in metallic properties, often referred to as "non- Fermi liquid behavior", that is 
induced over a wide range of temperatures above the quantum critical point (QCP), together with 
the marked tendency to develop superconductivity in the vicinity of such "quantum critical points" 
has given rise to a resurgence of interest in heavy fermion materials. 

The experimental realization of quantum criticality returns us to central questions left unan- 
swered since the first discovery of heavy fermion compounds. In particular: 

• What is the fate of the Landau quasiparticle when interactions become so large that the 
ground state is no longer adiabatically connected to a non-interacting system? 

• What is the mechanism by which the antiferromagnet transforms into the heavy electron 
state? Is there a break-down of the Kondo effect, revealing local moments at the quantum 
phase transition, or is the transition better regarded as a spin density wave transition? 



20031 ). a material with a 



Fig. [32j illustrates quantum criticality in YhRho Sio nCusteis et al\ . 
90mK magnetic transition that can be tuned continuously to zero by a modest magnetic field. 
In wedge-shaped regions either side of the transition, the resistivity displays the dependence 
p{T) = po + AT'^ (colorcoded blue) that is the hallmark of Fermi liquid behavior. Yet in a tornado 
shaped region that stretches far above the QCP to about 20K, the resistivity follows a linear 
dependence over more than three decades. The QCP thus represents a kind of "singularity" in the 
material phase diagram. 

Experimentally, quantum critical heavy electron materials fall between two extreme limits that 
I shall call "hard" and "soft" quantum criticality. "Soft" quantum critical systems are moderately 
well described in terms quasiparticles interacting wi th the soft quantum sp i n fluc tuations created 



by a spin-density wave instability. Theory predicts (|Moriva and Kawabata 



I973I I that in a three 



dimensional metal, the quantum spin density wave fluctuations give rise to a weak yT singularity 
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FIG. 32 (a) Color coded plot of the logarithmic derivative of resistivity d In p/ d In T after (jCusters et al 



20031 ) fb') Resistivity of YbRh2Si2 in zero magnetic field, after (jTrovarelli et a/.l . 120001 ). Inset shows logarith- 
mic derivative of resistivity 

in the low temperature behavior of the specific heat coefficient 

^ = 70 -liVf 



Examples of such behavior include CeNi2Ge2 (jGrosche et al 



200C 



Kiichler et al 



20031 ') chemi- 



cally doped Ce2-xLaxRu2Si2 and "A"-type antiferromagnetic phases of CeCu2Si2 at a pressure- 
tuned QCP. 

At the other extreme, in "hard" quantum critical heavy materials, many aspects of the 
physics appear consistent with a break-down of the Kondo effect associated with a re- 
localization of the f-electrons into ordered, ordered local moments beyond the QCP. Some of 
the most hea vily studied examples of this behavior occur in the chemically tuiied QC P in 



CeCuQ^x^Ux (jvon 



YbRh2Si2-xGex 
s of YbRh2Si2 



Niklowitz et al. 



Lohnevsenl. 



1996; 



Custers et al 



Trovarelli et al. 



von L5hnevsen et al. 



2003 



Gegenwart et al 



200fl) and YbAgGe{ 



1994 : 



Schroeder et al. 



1998 



2000). and 



20051'! an d the field-tuned tuned QCP 



Bud'ko et al 



20041 . 



2005 



Fak et al 



200,51 : 



20061 ) Hallmarks of hard quantum criticality include a logarithimically diverging 



specific heat coefficient at the QCP, 



and a quasi-linear resistivity 



To \t 



p{T) ~ r^+'' 



(136) 



(137) 



where t] is in the range — 0.2. The most impressive results to date have been observed at field- 
tuned quantum critical points in YbRh2Si2 and CeCoIn^, where linear resistivity has been seen 
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to extend over more t han two decades of t e mper a ture at the fie l d-tun ed quantum critical point 



(jPaglione et al. 



2001 



2006 



Ronning et al 



2006 



Steglich et al. 



19761 ). Over the range where 



linear , where the ratio between the change in the size of the resistivity Ap to the zero temperature 
(impurity driven) resistivity po 



CeCoIn^ is particular 



flow along the c-axis (jTanatar et al 



Ap/po » 1. 



y interesting, f or in this case, this resistance ratio exceeds 10^ for current 



20061 ). This observation excludes any explanation of that 



attributes t^ 



re unu s ual re sistivity to an interplay between spin-fluctuation scattering and impurity 



scattering (jRosch 



19991 ). Mysteriously, CeCoIn^ als o exhibits a T^^"^ re sistivity for resistivity 



for c urrent flow in the bas al plane below about 2K (jTanatar et al. 



Pines ([Nakasuii et al. 



2006 ). Nakasuji, Fisk and 



2004 ) have proposed that this kind of behavior may derive from a two fluid 
character to the underlying conduction fluid. 

In quantum critical YbRh2Si2~xGex, the specific heat coefficient develops a divergence 
at the very lowest temperature. In the approach to a quantum critical point, Fermi liquid behavior 
is confined to an ever-narrowing range of temperature. Moreover, both the linear co efficient of the 



speci f ic heat and the the qu adratic coefficient A of the resistivity appear to diverge (jEstrela et al 



2002 



Trovarelli et al 



2000). Taken together, these results suggests that the Fermi temperature 



renormalizes to zero and the quasiparticle effective masses diverge 



m 

0, — 

m 



oo 



(138) 



at the QCP of these three dimensional materials. A central property of the Landau quasiparti- 
cle, is the existence of a finite overlap "Z", or "wavefunction renormalization" between a single 
quasiparticle state, denoted by |qp~) and a bare electron state denoted by |e~) = c|^^|0), 

m. 

z 



{e |qp 



(139) 



If the quasiparticle mass diverges, the overlap between the quasiparticle and the electron state from 
which it is derived is driven to zero, signaling a comple te break-down in the quasiparticle concept 



at a "hard" quantum critical point (jVarma et al 



2002 ) 



Table. C. Selected Heavy Fermion compounds with quantum critical points. 
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Compound 


Xc/Hc 


Cv 

T 




r(r) = ^ 


Other 


Ref. 




CeCuQ-x^^x 


Xc = 0.1 




T + c 


- 




[1] 


"hard" 


YbRh2Si2 


5^11 = 0.66r 




T 


- 


Jump in Hall 
constant. 


[2] 




YbRh2Si2~xGex 


Xc = 0.1 




T 


2^-0.7 




[3] 




YbAgGe 


Bc\\ = 9-t 
Bc± = 5T 




T 


- 


NFL over 
range of fields 


[4] 




CeCoIn^ 


Bc = 5T 




1 2^1.5 




jjQ ± , 

Pab OC 


[5] 


"soft" 


CeNi2Ge2 


Pc = 


70 - 71 


2^1.2-1.5 


T-1 




[6] 





von Lohnevsen . 1 99t : 


von Lohnevsen e 


al. 


1994 


■ Schroeder et al. 


199^.200( 


),[2] 


Paschen et al. 


2004 


Trovarelli et al.. 


200C). [31 iCusters et al.. 


200c : Gesenwart et al. . i 


m5). [41 1 


Bud'ko et al. 2004. 


2005 


Fak et al.. 


2005. 


Niklowitz et al. 


200t ) [51 


Pat 


rlione et al.. 200c. 200t 


: Ronnins et al.. 


2006 : Tanatar et al. 


200( 


),m 


Grosche et al. 


200(.: 


Kiichler et al.. 


200c 


). 



Table C. shows a tabulation of selected quantum critical materials. One interesting variable 
that exhibits singular behavior at both hard and soft quantum critical points, is the Griineisen 
ratio. This quantity, defined as the ratio 

1 



a _ 1 91nr 
~C~~V ~dP~ 



OC 



of the thermal expansion coefficient a = to the specific heat C, diverges at a QCP. The 

Griineisen ratio is a sensitive measure of the rapid acquisition of entropy on warming away from 
QCP. Theory predicts that e = 1 at a 3D SDW QCP , as seen in CeNi2Ge2. In the "hard" 



quantum critical materi a.1 YhRh^Si^-x Ge 



density wave instability ([Kiichler et al. 



2003|). 



0.7 indicates a serious departure from a 3D spin 



B. Quantum, versus classical criticaltiy. 



Figure [551 illustrates some key distinctions between classical and quantum criticality (|Sachdevl . 



20071 ). Passage through a classical second-order phase transition is achieved by tuning the temper- 



ature. Near the transition, the imminent arrival of order is signaled by the growth of droplets of 
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nascent order whose typical size ^ diverges at the critical point. Inside each droplet, fluctuations 
of the order parameter exhibit a universal powerlaw dependence on distance 

1 



(V'(x)V'(o)) 



rd-2+r] ' 



(x « 0- 



(140) 



Critical matter "forgets" about its microscopic origins, for its thermodynamics, scaling laws and 
correlation exponents associated with critical matter are so robust and universal that they recur 
in such diverse contexts as the Curie point of iron or the critical point of water. At a conventional 
critical point order parameter fluctuations are "classical" , for the characteristic energy of the critical 
modes flu){qQ), evaluated at a wavevector qq ~ S,^^, innevitably drops below the thermal energy 
huj{qQ) « k^Tc as ^ ^ oo. 



19751 ) and Hertz (jHertd . Il976l ) recog 



In the seventies various authors, notably Young (jYouna . 
nized that when if the transition temperature of a continuous phase transition can be depressed 
to zero, the critical modes become quantum mechanical in nature. The partition function for a 



quantum phase transition is described by a F eynma.n integral over order 



{^(x,r)} in both space and imaginary time (|Hertzj . 



1973; 



Sachdevl . 



parameter configurations 



2003) 



Z. 



quantum 



(141) 



space-time configs 



where the action 



S[^^J] 



dT 



d'^xL[i;{x,T)], 



(142) 



contains an integral of the Lagrangian L over an infinite range in space, but a Enite time interval 

If 



Near a QCP, bubbles of quantum critical matter form within a metal, with finite size and 
duration (Fig. [33]) . These two quantities diverge as the quantum crit ical point is approached, 



but the rates of divergence are related by a dynamical critical exponent (jHertz 



19761 ). 



One of the consequences of this scaling behavior, is that time counts as z spatial dimensions, 
[t] = [L^] in general. 

At a classical critical point, temperature is a tuning parameter that takes one through the 
transition. The role of temperature is fundamentally different at a quantum critical point: for it 
sets the scale ~ 1/T in the time direction, introducing a finite size correction to the quantum 
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QUANTUM P 

FIG. 33 Contrasting classical and quantum criticality in heavy electron systems. At a QCP an external 
parameter P, such as pressure or magnetic field replaces temperature as the "tuning parameter" . Temper- 
ature assumes the new role of a finite size cutoff 1^ cx 1/T on the temporal extent of quantum fluctuations. 
(A) Quantum critical regime, where Ir < S,tau probes the interior of the quantum critical matter (B) Fermi 
liquid regime, where l^- > S,r, where like soda, bubbles of quantum critical matter fleetingly form within a 
Fermi liquid that is paramagnetic (Bl), or antiferromagnetically ordered (B2). 

critical point. When the temperature is raised, Ir reduces and the quantum fluctuations are probed 
on shorter and shorter time-scales. There are then two regimes to the phase diagram, 

(A) Quantum critical: Ir « (143) 
where the physics probes the "interior" of the quantum critical bubbles, and 

(B) Fermi liquid/AFM Ir » (144) 

where the physics probes the quantum fluid "outside" the quantum critical bubbles. The quan- 
tum fluid that forms in this region is a sort of "quantum soda", containing short-lived bubbles 
of quantum critical matter surrounded by a paramagnetic (Bl) or antiferromagnetically ordered 
(B2) Landau Fermi liquid. Unlike a classical phase transition, in which the critical fluctuations 
are confined to a narrow region either side of the transition, in a quantum critical region (A), 
fluctuations persist up to tempe ratures where Ir becomes com parable the with the microscopic 



short time cut-off in the problem (jKopp and Chakravartv 



20051 ) (which for heavy electron systems 



is most likely, the single-ion Kondo temperature ~ Ti/Tk)- 
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FIG. 34 (a) Hall crossover line for sudden evolution of Hall constant in YhRhiSi2 after (jPaschen tt al 



20041 ). (b) Sudden change in dHvA freque ncies and divergence of quasiparticle effective masses at pressure 



20051) . 



tuned, finite field QCP in CeRhln^ after (jShishido et al 
C. Signs of a new universality. 



The discovery of quantum criticality in heavy electron systems raises an alluring possibility of 
quantum critical matter, a universal state of matter that like its classical counterpart, forgets its 
microscopic, chemical and electronic origins. There are three pieces of evidence that are particularly 
fascinating in this respect: 

1. Scale invariance, as characterized b y E/T scaling in the quantum -critical inelastic spin fluc- 



tuations observed in CeCui^xAux ( Schroeder et al. 



19981 . 



2000). {x = Xc = 0.016), 



AoiE,T) = -F{E/T) 
where a ~ 0.75 and .Ffxl oc (1 — Similar behavior has also been seen in powder samples 



of UCu^-xPdx (jAronson et al. 



19951). 



2. A ju mp in the Hall consta nt of YbRh2Si2 when field tuned through its quantum critical 



point (jPaschen et al. 



2004 1. (see Fig. [31(a)). 



3. A sudden change in the area of the extremal Fermi surface orbits o bserved by de Haas van 

■ ■ 12005 ). (see Fig. [341(b)). 



Alphen at a pressure tuned QCP in CeRhin^ (|Shishido et al 



Features (2) and (3) suggest that the Fermi surface jumps from a "small" to "large" Fermi surface 
as the magnetic order is lost, as if the phase shift associated with the Kondo effect collapses to zero 
at the critical point, as if the f-component of the electron fluid Mott-localizes at the transition. 
To reconcile a sudden change in the Fermi surface with a second-order phase transition, we are 
actually forced to infer that the quasiparticle weights vanish at the QCP. 
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These features are quite incompatible with a spin density wave QCP. In a spin density wave 
scenario, the Fermi surface and Hah constant are expected to evolve continuously through a QCP. 
Moreover, in an SDW description, the dynamical critical exponent is z = 2 so time counts as z = 2 
dimensions in the scaling theory, and the effective dimensionality D^ff = d + 2 > 4 lies above the 
upper critical dimension, where mean-field theory is applicable and scale invariant behavior is no 
longer expected. 

These observations have ignited a ferment of theoretical interest in the nature of heavy fermion 
criticality. We conclude with a brief discussion of some of the competing ideas currently under 
consideration. 



1. Local quantum criticality 



One of the intriguing observations (jSchroeder et al.l . 119981 ) in CcCuq^xAux, is that the uniform 
magnetic susceptibility, ^ T°- + C, a = 0.75 displays the same powerlaw dependence on 
temperature observed in the inelastic neutron sca ttering at the critical wavevector Qq. A more 



detailed set of measurements by Schroeder et al. (jSchroeder et al. 



2000) revealed that the scale- 



invariant component of the dynamical spin susceptibility appears to be momentum independent, 



X 



-1 



[c^,E)=T'^[^{E/T)]+x,\c^) 



This behavior suggests that the critic al behavior associated with the heavy 



some kind of local critical excitation (|Colemanl . 



19991 : 



Schroeder et al. 



(145) 



"ermion QCP contains 



1998|) 



One possibili t y, is t hat this local cr itical excitation is the spin itself, so that (j Coleman 



Sachdev and Yd . 



1993 



Sengupta , 



2000) 



19991 : 



{S{r)S{T')) 



(146) 



(r-r')2-^' 

is a power-law, but where e 7^ signals non- Fermi liquid behavior. This is t he basis of the ' "local 
quantum crit i cality " theory developed by Si, Ingersent, Rabello and Smith (ISi et al. 



Smith and Sil . 



2001 



20031; 



2000). This theory requires that the local spin susceptibility xioc = Z]q x(ct) '^)aj=o 



diverge s at a heavy fermion QCP. Using an ext ension of the methods of dynamical mean-field 



theory ([Georges et al. 



1996 



Kotliar et al. 



20061 ). Si et al. find that it is possible to account for 
the local scaling form of the dynamical su sceptibility, obta i ning e xponents that are consistent with 



the observed properties of CeCu^-xAux (IGrempel and Sil . 



20031). 



However, there are some significant difficulties with this theory. First, as a local theory, the 
quantum critical fixed point of this model is expected to possess a finite zero-point entropy per spin. 
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a feature that is to date, inconsistent with thermodynamic measurements (jCusters et all 120031 ) . 
Second, the requirement of a divergence in the local spin susceptibility imposes the requirement 
that the surrounding spin fluid behaves as layers of decoupled two dimensional spin fluids. By 
expanding Xo ^(q) (|145p about the critical wavevector Q, one finds that the singular temperature 
dependence in the local susceptibility is given by 

XUT) ~ / - (147) 

requiring that d < 2. 

My sense, is that the validity of the original scaling by Schroeder et al still stands and that 
these difficulties stem from a mis-identification of the critical local modes driving the scaling seen 
by neutrons. One possibility, for example, is that the right soft variables are not spin per-se, but 
the fluctuations of the phase shift associated with the Kondo effect. This might open up the way 
to an alternative formulation of local criticality. 



2. Quasiparticle fractionalization and deconfined criticality. 



One of the competing sets of ideas under consideration at present, is the idea that in the 
process of localizing into an ordered magnetic moment, the composite heavy electron breaks up 
into constituent spin and charge components. Loosely speaking. 



(148) 



where represents a neutral spin-1/2 excitation, or "spinon" . This has led to proposals 



( Coleman et al 



2001 



Pepin 



20051 : 



Senthil et al. 



20031 ) that gapless spinous develop at the quan- 



tum critical point. This idea is faced with a conundrum, for even if free neutral spin-1/2 excitations 
can exist at the QCP, they must surely be confined as one tunes away from this point, back into the 



Fermi liquid. According to t^ 



le model of "decon fined criticality" proposed by Senthil, Vishwanath, 



Balent, Sachdev and Fisher (jSenthil et al 



20041 ) , the spinon confinement scale .^2 introduces a sec- 
ond diverging length scale to the phase transition, where ^2 diverges more rapidly to infinity than 
^1. One possible realization of this proposal is the quantum melting of two dimensional 5 = 1/2 
Heisenberg antiferromagnet, where the smaller correlation length ^1 is associated with the transi- 
tion from antiferromagnet to spin liquid, and the second correlation length ^2 is associated with 
the confinement of spinous to form a valence bond solid (Fig. I35|) . 

It is not yet clear how this scenario will play out for heavy electron systems. Senthil, Sachdev 



and Vojta (jSenthil et al 



2OO5I ) have proposed that in a heavy electron system, the intermediate 
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Valence bond 
Solid 



FIG. 35 "Deconfined criticality" (jSenthil et 



2004 ). The quantum critical droplet is defined by two 



divergent length scales - governing the spin correlation length, ^2 on which the spinons confine, in the 
case of the Heisenberg model, to form a valence bond solid. 

spin liquid state may involve a Fermi surface of neutral (fermionic) spinons co-existing with a 
small Fermi surface of conduction electrons which they call an FL* state. In this scenario, the 
quantum critical point involves an instability of the heavy electron fluid to the FL* state, which is 
subsequently unstable to a ntiferromagneti s m. Re cent work suggests that Hall constant can indeed 



jump at such a transition (jColeman et al 



2005a|). 



3. Schwinger Bosons 



A final approach to quantum criticality currently under development, attempts to forge a 
kind of "spherical model" for the antiferromagnetic quantum critical point through the use of 
a large expansion i r i whi c h the spin is described by S chwinger bosons, rather than fermions 



(jArovas and Auerbachl . 



1988 



Parcollet and Georgesl . 



1993), 



Sab = bah — ^ab^ 

where the spin S of the moment is determined by the constraint ni, = 2S on the total num- 
ber of bosons per site. Schw inger bosons are well-suited to describe low dimensional magnetism 



(jArovas and Auerbachl . 



19881 ). However, unlike fermions, only one boson can enter a Kondo singlet. 



To obtain an energy that grows with A^, Parcollet and Georges proposed a new class of large N 
expansion based around the multi-channel Kondo model with K channels (jParcollet and Georgesl . 
19971 ). where k = K/N is kept fixed. The Kondo interaction takes the form 



H. 



int 



Jk \ - 



u=l,K, a, P 



(149) 
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where the channel index v runs from one to K. When written in terms of Schwinger bosons, this 
interaction can be factorized in terms of a charged, but spinless exchange fermion Xv ("holon"), as 
follows 



Jk 



(150) 



arcollet and Georges origin ally used this method to study the over-screened Kondo model 



Parcollet and Georges 



1993), where K > 2S. 

Recently, it has proved possible to find the Fermi liquid large solutions to the fully screened 
Kondo im purity model, where the number of channe l s is co mmensurate with the number of bosons 



{K = 2S) (ILebanon and Colemanl . 



2007 



Rech et al. 



20061 ) . One of the intriguing features of these 



solutions, is the presence of a gap for spinon excitations, roughly comparable with the Kondo tem- 
perature. Once antiferromagnetic interactions are introduced, the spinous pair-condense, forming a 
state with a large Ferni i surface, but one that co-exists with gapped spinon (and holon) excitations 



(Coleman et al 



2005b|). 



The gauge symmetry associated with these particles guarantees that if the gap for the spinon 
goes to zero continuously, then the gap for the holon must also go to zero. This raises the possibility 
that gapless charge degrees of freedom may develop at the very same time as magnetism (Fig. [36] 
). In the two impurity model, Rech et al have recently shown that the large N solution contains 
a "Jones Varma" quantum critical point where a static valence bond forms between the Kondo 
impurities. At this poi nt, the holon and spin on excitations become gapless. Based on this result, 



Lebanon and Coleman (Lebanon et al. 



20061 ) have recently proposed that the holon spectrum may 
become gapless at the magnetic quantum critical point (Fig. [36l ) in three dimensions. 



VI. CONCLUSIONS AND OPEN QUESTIONS 

I shall end this chapter with a brief list of open questions in the theory of heavy fermions. 

1. To what extent does the mass enhancement in heavy electron materials owe its size to the 
vicinity to a nearby quantum phase transitions? 

2. What is the microscopic origin of heavy fermion superconductivity and the in extreme cases 
UBeis and PuCoGa^ how does the pairing relate to both spin quenching and the Kondo 
effect? 
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FIG. 36 Proposed phase diagram for the large N hmit of the two impurity and Kendo lattice models. 
Background- the two impurity model, showing contours of constant entropy as a function of tem perature and 
the ratio of the Kondo temperature to Heisenberg coupling constant, after (jRech et a/.l . 120061 ) . Foreground, 
proposed phase diagram of the fully screened, multi-channel Kondo lattice, where S is the spin of the 
impurity. At small S*, there is a phase transition between a spin liquid and heavy electron phase. At large 
5, a phase transition between the antiferromagnet and heavy electron phase. If this phase transition is 
continuous in the large N limit, th en both the spinon and holon gap are likely to close at the quantum 
critical point (jLebanon et aZ.1 . 120061 ) . 



3. What is the origin of the linear resistivity and the logarithmic divergence of the specific heat 
at a "hard" heavy electron quantum critical point? 

4. What happens to magnetic interactions in a Kondo insulator, and why do they appear to 
vanish? 



5. In what new ways can the physics of heavy electron systems be interfaced with the tremen- 



dous current developments in mesoscopics? 



feature of Coulomb blockaded quantum dots (jKouwenhoven and Glazman 



he Kondo effect is by now a well- established 



2001 



), but there 



may be many other ways in which we can learn about local moment physics from mesoscopic 
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experiments. Is is possible, for example, to observe voltage driven quantum phase transitions 
in a mesoscopic heavy electron wire? This is an area ripe with potential. 

It should be evident that I believe there is tremendous prospect for concrete progress on many of 
these issues in the near future. I hope that in some ways, I have whet your appetite enough to 
encourage you too try your hand at their at their future solution. 
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